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ABSTRACT 


A frequency dependent hybrid mode analysis of genexal 
coplanar transmission lines is presented. A Fourier trans- 
form ЕСЕ is applied and the resulting expressions are 
evaluated numerically using the method of moments. Shielded 
slot line resonators as well as dispersion characteristic and 
characteristic impedance of open and closed boundary slot 
lines, coupled slot lines and coplanar strip lines are obtained. 
An extension of this technique is used in a perturbational 
analysis of slot lines on a ferrite substrate. 

The theoretical performances of various slot line direc- 

· tional couplers are predicted. Numerical results are compared 
with experiments in various cases, where excellent agreement 


has been obtained. 
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1.0 INTRODUCTION 


The study of electromagnetic energy transmission is 
but one important area in microwave engineering, where the 
electromagnetic waves are travelling through some trans- 
mission medium, which provides the link between the trans- 
mitting and the receiving part of a transmission system. 
The transmission medium may be the free space or it may be 
a guidance structure which will direct the energy flow into 
a certain prescribed direction. This structure will be in 
general a boundary where the surrounding fields interface 
with the surface currents and charges on this boundary which 
will be in most cases an electric conductor, although in 
some cases a dielectric material may be used instead. Some 
standard examples for guided wave structures for microwave 
applications are the cpen wire transmission lines, coaxial 


cables and rectangular or circular waveguides. 


1.1] TRANSMISSION LINES FOR INTEGRATED CIRCUIT APPLICATIONS 

In most cases of designing a transmission link the deci- 
sion of which guidance structure should be used will be based 
upon various determining factors which are sought to be opti- 
mal in some sense. Some examples of these factors are the 
cost of design and maintenance, the transmitted power level 
and the frequency band of operation, to name only a few. 


While most of these examples given so far apply mainly to 
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large scale transmission SyStems, recent advances and 
increasing use of integrated circuit technology at microwave 
frequencies set different standards where new transmission 
lines are required that are planar in nature. Several rea- 
sons may have led to this ever increased interest in planar 
transmission lines; e.g. the use of microwave circuits in 
Satellite communications requires light weight and small 
Sizes; another reason is that the etching of a planar trans- 
mission line on a dielectric or ferrite substrate is just 
another production step within the fabrication of a wholly 
automatically fabricated integrated circuit module. The most 
commonly used transmission line to date in microwave IC's 
is the microstrip line shown in Figure 1.1/1. It was first 
proposed in 1955 by M. Arditi [Ref. 1] and has since then 
received considerable attention by innumerable investigators. 
Although the investigations of single and coupled micro- 
Strip lines continue to be the subject of more sophisticated 
research today only recently the interest shifted towards 
investigations of coplarar transmission lines. Some examples 
of these lines are shown in Figures 1.1/2-4. These lines 
offer the advantage of being more easily adaptable to shunt 
element connections without the need to penetrate the die- 
lectric substrate as in the case of microstrip lines. Active 
devices as well as passive elements may be fabricated in 
either series or shunt configuration directly on to these 
lines in either hybrid or monolithic fabrication techniques 


mer. 2]. | 
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Figure 1.1/4 Parallel Coplanar Strip Transmission Line 
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The slot line in Figure 1.1/2 consists of a narrow gap 
between two conducting surfaces which are on one side of a 
dielectric substrate. A fairly high permittivity of this 
substrate confines the waves travelling on this structure to 
the region of this slot such that the main mode of propaga- 
Gon is similar to a TE 


“10 


Figure 1.1/3 shows two adjacent slot lines parallel to 


rectangular waveguide mode [Ref. 3]. 


each other and of equal width, where the coupling effect 
between these lines may be advantageous in the design of 
microwave filters and couplers. ТЕ this configuration is 
used as a single transmission line such that both outer con- 
ductors are the ground electrodes, this line is called a 
coplanar waveguide (CPW) according to Wen [Ref. 4]. 

Both iines use a fairly broad area of conductor surface 
which may serve as the ground plane for a microstrip placed 
onto the opposite side of the substrate, thus coupling from 
one transmission line to the other can be easily achieved 
Mer. 5]. 

Another promising property of these two lines is the fact 
that both provide an area where the RF magnetic field is 
elliptically polarized which makes these lines suitable for 
nonreciprocal applications when the dielectric is replaced by 
осете substrate. This elliptic polarization is in con- 
trast to microstrip lines where a similar effect can be achieved 
only by using meander type lines [Ref. 2]. The coplanar strip 


line shown in Fig. 1.1/4 consists of two strip conductors on 





a dielectric substrate with no ground plane and can be 


thought of as a dual to the coplanar waveguide [Ref. 4]. 


1.2 MATHEMATICAL ANALYSES OF COPLANAR TRANSMISSION LINES 

In the desiyn of a transmission line the two important 
line parameters are the characteristic impedance and the 
propagation or phase velocity. 

A standard approach in the analysis of these lines is to 
consider the line as having distributed parameters of capaci- 
tance ard inductance and to solve for these in terms of the 
physical dimensions. This solution assumes a TEM mode of 
propagation Since static fields are assumed for the calcula- 
tion of the line reactances. However it can be easily proved 
that the inhomogeneous guidance structures shown in Figures 
1.1/1-4 strictly cannot support a TEM mode. Assume waves 
are propagating along a non homogeneous guiding system such 
as a microstrip or a slot line with a propagating term 
е (94-2) where w is the radian frequency, t the time and 
У = о + jR the propagation constant assumed from now on to 
be pure imaginary (lossiess case, a=0), ј = m 

The transverse field components can be related to the 
longitudinal field components through an algebraic manipula- 


tion of Maxwell's curl equations to yield 


H cuml 8 (we р № diu 
xi Ze | | оу "3x 


к. ~B 
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If the TEM mode is assumed E, and H, are non existing 
throughout the guidance structure. A solution other than 
the trivial solution clearly can only exist if k; -B = 0. 
Since for the inhomogeneous guidance system ß is the same 
Sr all regions but for certain 173, E the TEM mode 
condition cannot be met. 

But although a planar transmission line such as the 
microstrip line has an inhomogeneous cross-section a TEM 
or often called "quasi-static" mode approximation can be used 
Since most of the fields are confined to the dielectric be- 
tween the two conductors. Recent investigations show however 
that this is only valid for frequencies up to approximately 
5 GHz [Ref. 6]. For coplanar transmission lines a quasi- 
Static approximation seems no longer suitable for microwave 
frequencies in the GHz region since the waves are not as 
closely bound in the dielectric. Investigations of coplanar 
transmission lines thus should start from the general field 


theory to include the hybrid mode which provides a frequency 
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dependent analysis, This means that in general a solution 
to Helmholtz's equation rather than Poisson's equation is 
to be found with inhomogeneous boundary conditions due to 
the interaction between fields and guiding structure. This 
is not an easy task in a mathematical sense. A Green's 
function solution e.g. should be feasible at least theore- 
tically but has not been reported in the literature so far, 
probably due to the immense complexity of the problem. 

In his analysis of the slot line Cohn [Ref. 3] introduced 
electric and magentic boundary walls which permitted the 
Slot line configuration to be treated aS a rectangular wave- 
guide problem. The set of er waveguide modes was 
Е саз а set Of Orthogonal functions; in the final analysis 
the fields were approximated by the ТЕ 10 rectangular waveguide 
mode. 

Wen [Ref. 4] investigated the coplanar waveguide and the 
coplanar strips with the assumption that the dielectric 
substrate is thick enough to be considered infinite for 
conformal mapping purposes. For large values of the permit- 
tivity this may be a practical assumption but it appears to 
be impractical for relatively small values or for thin 
substrates. An alteration of this method used by Wen is 
given by Davis et al. [Ref. 7] and takes the finite thickness 
of the substrate into account. The coplanar strip line, 
however, is not treated. Furthermore, both methods lack any 
frequency dependent information on the behavior of the line 


parameters since a quasi-static approximation is used. 
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Recently, Itoh and Mittra [Ref. 8] proposed a new method 
for calculating the dispersion characteristic of a single 
slot line. Here Helmholtz's equations for the scalar 
potential functions [Ref. 9] are transformed into the Fourier 
domain and are solved for in accordance with the given 
boundary conditions. The method of moments is then used to 
solve for the propagation constant. This approach showed 
good agreement with results from Cohn [Ref. 3]. However no 
analysis of the characteristic impedance has been reported 
in Reference 8. 

This spectral domain approach was used by Luna [Ref. 10] 
to compute the dispersion characteristic of coplanar strips. 
No further analysis of this transmission line has appeared 


еспе literature. 


1.3 UNSOLVED COPLANAR TRANSMISSION LINE PROBLEMS 

Besides the known parameters of coplanar transmission 
lines given in Section 1.2 there remain several problems open 
for analysis which will be shortly listed here and are the 
subject of this thesis. Although the slot line was studied 
in 1968 by Cohn [Ref. 3] with a frequency dependent analysis 
and the coplanar waveguide structure was analyzed by Wen 
and Davis et al. [Refs. 4,7] by a quasi-static method, no 
complete analysis of two coupled slot lines has appeared yet, 
although originally planned by Cohn. Line parameters which 
describe the coupling effect between these lines are important 


in the design of directional couplers end in filter design. 
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Moreover, in the design of slot line filters, the application 
of resonant slot structures may be desirable and these have 
not yet been studied analytically at all. 

Another unsolved problem is the shielded slot line 
Structure. Since waves traveling on a slot line are not 
closely bound to the dielectric next to the slot Bi extend 
somewhat into the air regions it seems important to study 
the behavior of the line parameters under the influence of 
nearby obstacles such as an enclosing or shielding structure. 
These results could be of some significance in the design of 
microwave IC modules. 

A further investigation is needed for a slot line on a 
ferrite substrate. Although recently a publication by 
Kitasawa 2. al. [Ref. 11] discussed this problem, the analysis 
presented is restricted since the operating frequency is 
assumed to be far a the resonance frequency of the 
ferrite. No analysis about attenuation of the waves on the 
ferrite is mentioned. 

Finally, although the frequency dependent dispersion 
characteristics of coplanar strips were analyzed by Luna 
[Ref. 10] a frequency dependent formulation of the charac- 


teristic impedance is still missing. 


1.4 CONTENTS AND SUMMARY 
Inhomogeneous waveguide structures such as coplanar 
transmission lines will not Support a TEM-mode of propagation 


at high frequencies. Therefore any analysis will assume a 





full or hybrid mode which may be thought of as being a 
superposition of TE and TM modes. This formulation will 
yield a frequency dependent phase velocity and characteristic 
impedance. 

This thesis is devoted to a hybrid mode analysis of 
coplanar transmission lines. In Section 2.0 open boundary 
structures are analyzed. Following the proposed method of 
Itoh and Mittra [Ref. 8] the dispersion characteristic of 
a Single slot line is formulated. The investigations are 
extended by comparing various possible approximations. Using 
Parseval's relation it is shown that the power flow on this 
Waveguide Structure can be computed in the spectral domain. 

The reSultant calculations for the characteristic 
impedance A compared with results from Reference 12. This 
is the first verification of Cohn's method by an entirely 
different numerical analysis and shows very close agreement. 
The theory is extended for an analysis of two coupled slot 
lines of equal width, the results are new and very important 
for the microwave industry in the design of slot line couplers 
and. filters. A frequency dependent analysis of coplanar 
Waveguide is inciuded. 

An easy extension of the theory is then used to analyze 
the line parameters of a coplanar strip transmission line. 
The final part of Section 2.0shows a perturbational expression 
for the propagation constant of slot lines on a ferrite 


substrate. 
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In Section 3.0 various closed boundary structures are 
analyzed. The shielded slot transmission line is compared 
with the open boundary structure.  Resonating slot configura- 
tions are investigated such as a single resonant slot and 
two end coupled resonating slots. A byproduct of the 
resonant slot investigation is the end effect analysis of 
Slot lines. Theoretical results are compared vith experi- 
mental data from open boundary experiments and show good 
ERureement. 

Section 4.0 is devoted to an experimental verification 
of the numerical data for coupled slot line parameters. 
Using basic network theory the theoretical response of a 
slot line directional coupler can be computed. Due to the 
non-equal phase velocities in the two natural modes co- or 
contra-directional coupling can be achieved. Two slot line 
КО ега, a co- and a contra-directional, were built and 
their performance measured. 

In Section 5.0 remarks are made about the numerical analvsis 
and the computer programs. It is believed that results of 
this thesis are of interest to the microwave industry. 

Hence two of the most important computer programs used in 
this work are attached to the thesis and their use explained 
in this section. 

Section 6.0 shortly outlines an attempt to analyze a 
circular resonating aperture on a dielectric substrate and 
shows the difficulties encountered there. Suggestions for 


expansions of the presented work conclude this section. 
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2.0 OPEN BOUNDARY STRUCTURES 


In this section open boundary coplanar tranSmission lines 
are analysed where it is assumed that the lines extend infi- 
mitely in the z-direction in an infinite x-z plane. It is 
further assumed that conductors are negligibly thick ee 
dielectric and conductor are perfect. 

The electromagnetic fields surrounding the guidance struc- 
ture are derived in a most general way from hybrid modes where 
use is made cf the fact that any hybrid mode can be thought 
of being a superposition of TE and TM modes by the orthogon- 
ality properties of open boundary waveguide modes, proved by 
Collins [Ref. 13]. The derivation of ail equations is general 
and not а to any specific mode of propagation. Approxi- 
mations for the lowest order mode are made in the very end 
and can be easily modified to yield results for higher order 
modes. 

It is further assumed that the operating frequency is kept 
below the frequency for the onset of surface wave modes which 
could propagate in the dielectric-air medium for y « D (see 
Figure 1.1/2-4). 

This onset of surface wave modes for the slot line con- 
figuration e.g. can approximately be determined. Consider the 
region y « D and assume that the dielectric has the permit- 
LEY E) while for y < O the permittivity is that of free 


Space E_. All further assumptions stated above hold and the 
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existence of a slot in the metal is irrelevant.  Attwood 
in Reference 14 has derived the determinantal equation for 


TM Surface wave modes which yield the lowest order mode. 


€ 
This equation is given as „ tan(k_ „D) = - = k E where 
E B E 2 2 
= к Ва апа ken ar Ко ; Ko оС апа 
39 =: w enuo A bound for the lowest order cutoff wave- 


length can be found by approximating this transcendental 


equation. For К 2 tan(k D) < 0, Сап (Кр) 9 0 (246 





: T _ 227272 | 

Ko)D > 5 · But к» = к. В so that the following 

inequality holds k AD » Vk, *- 8? р > 5 . ТЕ do ema Sill 
Еи © 


the bound for the onset of surface waves is 
Ae db.) , 


2.l SINGLE SLOT-LINE 

The schematic of a single slot line and its orientation 
with respect to a three-dimensional rectangular coordinate 
system is shown in Figure 1.1/2. The substrate is assumed 
to be dielectric. Let this transmission line be divided into 


three regions in space and define these regions as 


region 1 FOr D- y 
2 ОЛ КО SUD 
3 ү < 0 
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Since the lossless case is essumed, the propagation constant 
Y- oat j8 is pure imaginary with o= 0. The following 
derivation, however, is perfectly oeneral and the above 


simplification is made onlv at the end. 


К У оп Ее ИЕ СЕ 
In the construction of solutions to Maxwell's field 
equations for a transmission line problem it is sufficient 
to assume two scalar potential functions, a scalar electric 
potential ¢ and a scalar magnetic potential y which may be 
both related to TM and TE modes respectively [Refs. 9,13]. 


The axial components of TM and TE modes are then 


а. 2 tyz _ 
Е, 7 Kei ф(х,ује (2-1) 
= 2 tyz _ 
and H, - Koi v(x,y)e (2-2) 
2 2 2 A CHR t % % » 
where к = к. КИЕ, і-12,3 Xdefimngcthesspatial region 


of the transmission line as stated in Section 2.1. 
Through Maxwell's curl equations the transverse field 
components are then determined by these axial components and 


can be given as 





36. эф. 
= M E MEC M 
E di Ser ШЫ y ) е | | (2250) 
36. д 
ы ai й Е уг 
• s ee + 1 — 
Ei (у 57 jeu x ) е (2-4) 
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ду. ab. 





м am - en yZ E 
Hoi Y IT + Jue, y ) e (2525) 
дї. Эф. 
= - oL — і yZ = 
На (y = JU£: 3c ) e (2-6) 


15] 72913 
where propagation in the z direction is assumed and the 
dependence of the scalar potentials on the variables x and 
y is understood. 
The scalar potential functions satisfy Helmholtz's 


equations in the three spatial regions, thus 


(V + К 


ху ci к о | (2-7) 


and D + Koj ) Y, = 0 (2-8) 
where E denotes the two-dimensional Laplacian operator 

in the transverse (x,y) direction. So equations (2-7) and 
(2-8) are second order partial differential equations (P.D.E.). 
A common practice in solving these P.D.E. which essentially 
present a boundary value problem is to transform these into 
Smainary differential equations (O.D.E.) by an integral 
transform technique [e.g. Ref. 15] which may then be solved 
analytically. The boundary alone гот (2-7) and (2-8) 

are intuitively obvious, Ж ее Ф; (х,у) + 0 for 

х,у > to and furthermore f lò; (x,y) | ax is bounded for 
finite energy reasons — me same holds if $6 is replaced by v — 


since the time average power flow is related to the scalar 


S 





potentials as shown in Reference 13 and is assumed to be 
finite. Although no rigorous proof will be given it is 
therefore assumed that the scalar potential functions can 
be transformed into the Fourier domain via a Fourier 
transformation defined as 


co 


Р(ф(х,у)) = Ф(а,у) = / ф(х,у)е2% 


ах (2--9) 


Tf this transform is applied to equations (2-7) and (2-8) 


the following O.D.E. are obtained: 


2 
д 2 
ду 
2 | 
д 2 2 
me te) = у. 1, (ау) (2-11) 
ду 
мһеге ү; = a* E ү? = kn 


Again for energy reasons certain solutions to these 
O.D.E. may be eliminated by boundary conditions in the three 
spatial БЕШІШӘШЕ СаСО әмбе y 5) (х,у) 0 апа во tla,y) +0. 


Тһе solutions to (2-10) and (2-11) may then be given 


as 
есу) 
фу (аљу) = Ас(оје 7 
$ (а,У) = B* (a) sinh YoY a C* (a) cosh Y Y (2-12) 
YY 
$a (ay) =D (a)e ` 
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and "Y, (Y-D) 


Y, (0,y) = А? (а)е 


$2 (а,У) B? (a) sinh 4 ap Ban Y Y (2-13) 


Y4 (o, y) Ay e 

The following observation about the solutions in 
region 2 must he emphasized. Any wave on this inhomogeneous 
wavecuide structure will partly travel through air and partly 
through the dielectric slab. It is then obvious that 
ENS B < k in the lossless case where Ку = к; = К 


О 2 
2-42 


Since -% <a <% , Ғор |а| < к. 8 


E 
then s < 0 
For these values of the transform variable a the solutions 
in (2-12) and (2-13) for region 2. will change from hyperbolic 
59 ч я 5 = . " " ЕЕ 2E 2 
Boctrigonometric functions since jos NEL (oS NK o OF 5 
The eight unknown coefficients A* (a) pim рћ (а) 
are not independent but can be related to each other through 
the continuity conditions of the field components at the 


interfaces between the three spatial regions. These continuity 


conditions in the spatial domain are for 


а= 0 





OA = в (5,0, 2) 


(2-14) 
Поа) = ЧК ое 
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ama. tor y = D 


E, 4 (x,D,2) = Е,2(х,р,2) 
в (де Ix| < 2 
E, y (X1D,2) = 
0 elsewhere 
(2-15) 
E, (x)e'* |х| < Е 
Еу (х,р,2) = 
0 elsewhere 
3,00е7% |х| > 7 
Hu c D, Z) = Н, (oc) = 
и 0 elsewhere 
jQ00eY* — |х| > 9 
Ну (х,р,2) e Ho (XDZ) = | 
0 elsewherc 


— 


The subscript t is an abbreviation for tangential 
x and z components to avoid writing essentially identical 
equations twice. Ez, z) and ы are the electric field 
components and the current densities in the x-z plane at 
E D in region l. 

Denoting the Fourier transforms of the x- and z- 


directed current density and electric field components by 


E (a) = FIRE, (х) ) Е. (а) = Р(Е,(х)} 


J.(0) 9 FOG, G0) 7 „(а) = F(G, G0) 
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and applying these transforms and the transforms of the 
continuity conditions to equations (2-12) and (2-13) an 
elimination of the unknown coefficients AC through pe by 
substitution is possible. The details of this lengthy 
algebraic manipulation are given in Appendix A. The results 


yield a set of coupled equations of the form 


Ny (0,8) М2 (а,В) E, (a) J, (a) 
= (2-16) 
М. (а,В) МА (а, В) E, (a) J, (a) 


where the elements of the N matrix can be thought of being 
the ЕЕ transforms of dyadic Green's function components 
[Ref. 13] by an application of the convolution integral in 
the space domain. The lossless case ү = jB is assumed in 
(2-16). 

Nothing seems to be gained so far since neither the 
electric field and current density components nor their 
transforms are known. The coupled equations in (2-16) can 
however be further simplified by the methods of moments 
(Ref. 16], applied in the spectral domain. Define the scalar 
product on the space of complex-valued functions of the real 
variable a over the domain -® < а < ® according to 
Reference 15 as 


со 


о = en (о) а (ва) ас 


ae OD 


35 





and choose as a suitable set of weighting functions the 
Fourier transforms of the electric field components, which 
means Galerkin's method [Refs. 9,16j. The moment method 


applied to equation (2-16) leads to 


[ 
© 


<Му{@,8)Ё (а),Ё (а)> + <N,(a,Bß)E_(a),E,(a)> 


(2-17) 


i 
o 


«N4(a,8) E. (a) ,E, (a)? + «N,(2,8)E, (a) ,E, (a)? 


That the right hand side is zero can be easilv verified by 
Parseval's relation [Ref. 15], observing that at y = D the 
presence of an electric field component excludes the presence 
Of a current density component and vice versa, thus electric 
field and current density are orthogonal to each other. 

Up to this point the formulation of the problem is 
exact Since no approximation has been made yet. Note that 
the electric field components and their transform although 
unknown can be expanded in a Set of complete kasis functions 


so that in general 
E, (x) = à a. e... (x) 
(218) 


and Е 09 = 1 b e (x) 


The following analysis is restricted to a one term approxima- 


Born for each of the series in (2-18) since it vas found that 
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Suficiente accuracy could be obtained for all practical 
purposes. So n = l in equation (2-18) and this subscript 
will be omitted from here on. 

With E (о) ~ Fle, (x) } and Е (о) ~ ге о) 


(2-17) then becomes 


а 7 мусба, В) [Е (а) |80 +ь Л Мр(а,В)Е (аЕ %(а)да = 0 
| (2-19) 
= 0 


ars 


м, (9,8) Е, (о)Е, *(а)да+ь Л м, (0,8) [Е (о) | aa 
со = OO 

The dispersion characteristic of the slot line can now be 

found by varying B such that the determinant of the coefficient 
matrix of (2-19) is zero for a given set of physical parameters 
at a desired frequency @ operation. 

The remaining question is what kind of basis functions 
to choose. The choice of this complete set of basis functions 
is arbitrary in a mathematical sense, since as long as this 
set is complete, any closed form of the fielä components can 
be represented by it. However, the rate of convergence of 
the series representation will depend on how well the first 
few terms approximate the closed form. This in general 
requires some a priori knowledge of the true distribution. 

In order to determine the sensitivity of this analysis 
to the choice of basis functions, an investigation of various 


One term approximations was made. One choice for the electric 
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field components which requires almost no insight into the 
field problem is of the form 


1 |х| < 7 
(1) 
E (x) = | (2520232) 
0 elsewhere 
W 
-1 - iS х < 0 
ES nu = 1 O <x < EH (2-20b) 
Z 2 
0 elsewhere 


fre FOurler transforms of these basis functions are then 


E CP (a) = = віп (2%) (2~21a) 
dl (а) = j 2 [1 ~ cos ($5 ] (2-21b) 
Another 


choice which certainly approximates the 
field more closely is 





(2) == lx] <> 
ex^ (x) = Vas? ке” @=22а) 
0 elsewnere 
W 
X (53 * — | x << 
ef) (x) (2-22b) 
0 elsewnere 





Ehe reason that (2-22) are Superior basis functions to 

(2-20) is that (2-22b} approaches zero smoothly as 

lx] >> апа |х| +0 . Furthermore (2-22a) appreximates the 

fields more closely near the edge of the slot than (2-20a). 
The Fourier transform of (2-22) can now be found in 


ENS tandard table of Fourier transforms (e.g. Ref. 17], 


(2) = aW 
E. (a) = m Bg (52 (2-23а) 
(2) S: V2 W oW Е 


where Bg and B, are Besse lii unctions of tne=mirst kind of 
order zero and two respectively. 

Call the use of x and ¿directed field components 
the second order approximation while the use cf an x directed 
component only is called a first crder approximation. Then 
a reduction in computation time is possible using this first 
order approximation since only опе integral in (2-18) has 
to be evaluated instead of four for the second order 
approximation. These two formulations of the dispersion 
problem were programmed into a digital computer language and 
the root 8 which rendered the determinant of the coefficient 
matrix to zero was found numerically by an iteration scheme. 
Mails o е толаша are given in Section 5. First and 
second order approximations using equations (2-21) and (2-23) 
are compared in Figure 2.1/1 with results from Cohn's method 


in Reference 12. Although no indication abovt the rate of 
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Figure 2.1/1 Comparison of Approximations for Dispersion 
Еее ОЕ а Single Slot Line 
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convergence and hence the accuracy of the two series for a 
general problem can be obtained, a comparison for this 
particular example shows that the largest deviation between 
the different approximations is less than 4%. As already 
noted above, the basis set (2-22) is vastly Superior to 
basis set (2-20) in a one term approximation for the field 
components due to the physically impossible discontinuities 
of equation (2-20b). 

A numerical investigation with regard to the magnitude 
of the x and z components of the elctric field was mađe, using 
the field equations (2-23) in the second order approximation 
for an arbitrary choice of W/D = 0.2 and Comm 20. Up Бо а 
value of normalized frequency D/A = 0.03 the x directed 
electric field component vas found to have a magnitude greater 
than ten times that of the z directed component. This provides 
justification for the use of the more efficient first order 
approximation. As one might expect, the magnitude of the 
z component of the electric field increases for larger values 
of normalized frequency. However, the agreement between this 
and Cohn's [Refs. 3,12] method is excellent for a first order 
approximation. The theory presented here is a verification 
of Cohn's theory and in a sense a check on the numerical 
EE Sic by Itoh and Mittra [Ref. 8] since it was derived 


independently with only their very brief letter for guidance. 
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212 sG@haracteristic Impedance 

The definition of the characteristic impedance for 
an ideal transmission line is uniquely given by static 
quantities because these lines support a TEM mode wave. On 
the other hand, in investigations of pipe-like waveguides 
the wave impedance is widely used [Ref. 13], however the 
impedances are different for TM and TE waves. Since the 
Slot line is neither a TEM transmission line nor does it 
support a pure TE or TM waveguide mode but rather supports 
hybrid modes, no unique definition of the characteristic 
impedance can be found. 

Since the choice of definition is arbitrary Cohn's 
very practical definition in Reference 3 is useä in the 
following derivation of the characteristic impedance for a 


slot line. Define 
7 = = = Оша | (2-24) 


where A is the slot voltage and Ра is the time average 


vg 
power flow on the slot line which according to Reference 13 


may be computed as 
р = Ф ве // Ё, хН,* a, ахау) (2—25) 
5 2 
Observe that although the hybrid mode can be fully described 


as a superposition of TE and TM modes, the evaluation cf the 


time average power flow is not just a superposition of the 
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single time average power flow carried by each single mode. 
A simple manipulation unveils that in addition to the single 
TE and TM power terms there are cross terms as shown below 
in (2-26). Assume ү = 38 , which means for z dependence 
mz a wave traveling in the -z direction. 

The transverse field components of the integrand 
in (2-25) can be expressed in terms of the scalar potential 
Rimnctions as shown in equations (2-3) through (2-6). Equation 


(2-25) becomes with this substitution and dependence cf the 


Scalar potentials cn x and y understcod 


poo -in / / (wer? + GR] * ocB LG + 


2 


(02 +7) (56 50 - 5$ 25)) ау аж (2-26) 
The limits of integration are infinite since the slot line 
is an open boundary structure, which makes the use of 
Parseval's reiation [Ref. 15] feasible. By this method 
(2-26) is transformed into the spectral domain where again a 
double integral is obtained with the variables of integration 


а апа у instead of x and y. 


со 


i 1] 9%, 2 906,2 
"ae p 4. Re E E [- о “Bwejo]* - ES - бес |с) 
| - jag? (o U 2% у 
o^Buu|v|^ = 3a8^to adus Y) 
я * 
QE Co AC s YI do dy (2-27) 
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where the dependence of the scalar potentials on variables 

a and y is implicitly understood. Besides any dependence 

on physical parameters and the frequency of operations, these 
Scalar potentials are furthermore dependent on the propagation 
Еопзгапс B, which is known from Section 2.1.1, and on the 
assumed field distributions as given by equations (2-12), 
(2-13) and the detailed derivation in Appendix A. Note that 
the functional dependence on the variable y may be removed 
by analytical integration,which changes the double integral 
to a single integral. For a numerical integration on a 
digital computer the integral (2-27) has to be evaluated 
separately for the three spatial regions since according to 
equations (2-12) and (2-13) the integrand assumes different 
forms Sec. The integrands for the three regions are given 
in Appendix B in equations (B-1) through (B-4). 

It should be noted that the amount of algebraic 
complexity in these equations although still Bom large 
could be decreased for numerical purposes by neglecting the 
z-directed electric field component which means that the 
Spatial phase shift between E, and Е, required no attention 
in this algebraic derivation. It is shown in the comparison 
peu Cohn's method [Refs. 3,12] that this use of a first 
order solution results in good agreement. 

A computer program was developed which first computes 
the dispersion characteristic and then the average power flow 
where the electxic field was approximated by (2-22a). To 


SOlve for the characteristic impedance in equation (2-24), 
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the slot voltage has to be computed. Since the slot voltage 


W 
2 
ШЕ = / E, (x) dx , this integration may be accomplished 
EN 
2 
analytically; depending on the assumed electric field 
approximations vo for (2-22a) for instance. A comparison 


of results obtained by this formulation of the characteristic 
impedance with values obtained by Cohn's method in Reference 
12 is given in Figure 2.1/2. The very close agreement for 
the two arbitrarily chosen sets of parameters substantiates 
the accuracy of both methods and is the first known verifica- 
tion of Cohn's method. 

For comparative purposes it is noted that the field 
distribution (2-20a) resulted in lower values for the 
characteristic impedance by not more than 10% from those 
obtained from (2-223). 

A byproduct of this analysis is a consideration of 
ee power distribution. The spatial division of the slot 
line in this numerical formulation allows an analysis of the 
amount of power floving in each of these three regions. In 
Figure 2.1/3 the ratio of the power in each recion over the 
total power carried on the structure is shown as a function 
of normalized frequency for two different permittivities. 
several facts can be learned from these graphs: 

a) The amount of power carried in each region is frequency 


dependent. As the frequency increases more power is 
састава through region 2, the dielectric region. 


b) Only a small amount of the total power is found in 


regon 3% 


с) More power is carried through the air regions if the 
permi.ttivity is lowered. 
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These results mav have been known previously but for the 
first time are here quantitatively analyzed by the means of 
the spectral domain formulation. 

The following investigation indicates, however, 
that a slight error due to the approximation in the first 
order exists. This example serves more illustrative purposes 
but might perhaps stimulate further research. With the above 
formulation of the time average power flow a graphical display 
of the power density for the slot line may be obtained. For 
an arbitrarily chosen set of parameters the power density 
in the spectral domain is shown in Figure 2.1/4a versus the 
normalized transform variable aD for the spatial region 1 
ле the substrate. The unit on the ordinate is omitted 
since it is insignificant in this illustration. 

From Cohn's zero order analysis in Reference 3 it is 
obvious that the power density is the largest close to the 
slot and it is also obvious that it does not change its Sign. 
TO verify this the Fourier inversion of the spectral domain 
power density has to be taken to obtain the power density in 
the space domain. Since the analysis iS numerical in nature 
inversion by a discrete Fourier transform technique (DFT) on 
a digital computer seemed appropriate. The result is 
displayed in Figure 2.1/4b. Clearly the power is the largest 
close to the slot and decays rapidly but contrary to any 
expectation changes its sign. It can be easily verified that 


this behavior is due to the small dip at aD - 0 for the power 








20 y 2.0 ” 
Сүр 





Еј 





A pet ge 


-10.0 -5.0 | 5 — 10.0 


X 
D 


Figure 2.1/4 Power Density Distribution in a) Spectral, 
b) Space Domain for Region 1 of a Single Slot 
Line with Bem 20, W/D = 1.0 and D/A = 0.02 
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density distribution in the spectral domain. This dip is 
believed to originate from the first order approximation of 
the electric field. Although a similar phenomenon can be 
observed in an investigation of the power densities in regions 


2 and 3 it is seen that the error introduced is quite small. 


EN COUPLED SLOT LINES 

With two transmission lines placed adjacent to each 
other such that electromagnetic coupling effects the fields 
on both line structures, the design of directional couplers 
and filters is possible. Figure 1.1/2 shows two parallel 
slot lines whose analysis is the subject of this section. 
There are many ways to describe tho coupling effect between 
two transmission lines. According to Reference 18, which is 
a compilation of some twenty-nine papers on coupled trans- 
mission lines, the most common formulation is to define the 
natural modes of the guided wave structure. Since in this 
analysis two siots of equal width and constant separation are 
considered, define the two natural modes as the even and odd 


mode such that 


Е. (х) EZER) for the even 


and E, (x) -E (7x) for the odd mode. 


The existing theory developed so far in Sections 2,1.1 and 
2.1.2 can now be easily extended for the analysis of line 


parameters for these two modes of the coupled slot lines. 
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It is observed that the dyadic Green's function transforms 
in equation (2-16) are functions of a, 6 , Є. р ала the 
operating frequency but are independent of the siot configura- 
tion. The physical dimensions of the slot enter into the 
calculation only through the assumed field distributions or 
basis functions. Thus, it is only necessary to modify 
equations (2-20) and (2-22) such that the mathematical formula- 
tion of the basis functions corresponds to the physicai 
configuration and field distributions of the coupled slots. 


AS an example Figure 2.2/1 shows the assumed eiectric field 


distribution 
1 > <х < S +W 
W2 S+W, 2 
(23) =- (x- ——) 
(2) _ 2 2 
e (x) = (2-28) 
x 
51 as 9 > 
м. 2 SD = 


0 elsewhere 


Ше modification for the z-directed electric field 
component is similar. The following short manipulation leads 
Кле Fourier transform of this distribution: 


The shifting property of Fourier transforms states that if 


Fí$(x)) » lo) = Л T dx then a I = „29Х0 5 (0) 


=» ОО 


ЕД 





x) 


even mode 
— -— Odd mode 








Figure 2.2/1 Assumed Electric Field Component in x-direction 
VENTO RE ven and Odd Mode on Coupled Slot Lines 
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The Fourier transform of (2-22a) is given in (2-23a) as 
EC (а). Шеп the Fourier transforms оғ the electric fiela 


in (2-28) become for the even mode 


ES SA 
aes ^ e. ^ae la) 
= 2 cos(a 523) . E ies (2-28a) 


and for the odd mode 


S+W 


(2) 
Е So 


er (a) = +j2 sin(a 


МА (2-28b) 
Another change to adapt the formulation of Section 2.1 

for the coupled line analysis is necessary in oquation (2-24); 

since the total time average power surrounding the transmission 


lines is now due to two lines 


я 
ауа о,е 
for the odd and even mode respectively. 

Beyond these no further major modifications are necessary 
to use the previously developed theory and conputer prograns. 
In the following discussicn only results using a field 
approximation (2-28a) and (2-28b) are presented, There are 


two reasons for this: a) the approximation of equations 


E22) was found to be superior to (2-20) in Section 2.1 
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b) the use of a second order approximation became prohibitive 
for the dispersion characteristic of coupled slots since the 
amount of computer time increased tremendously due to the 
necessity of using a more sophisticated numerical integration 
scheme (details in Section 5). Besides, the time average 
power computation uses only a first order approximation 
anyway. 

Figures 2.2/2a-b show ncw an example of the numerical 
results for the dispersion characteristics and the charac- 
teristic impedances, These results are basically new in the 
literature and only a partial comparison with other results 
is possible. A quantitative comparision e.g. between these 
odd mode line parameters and the quasi-static parameters 
Calculated by Davis et al. in Reference 7 was made and is 
shown in Figure 2.2/3 for W/D = 1.9 — or in the notation of 
Reference 7, t = S — at frequencies of 1, 3 and 5 GHz. Good 
agreement is found for the lower frequency range. Note that 
the characteristic impedance of a CPW is one half of Zoo (оаа 
mode characteristic impedance) for equal slot widths. 

A qualitative check on these results can be made by 
investigating the even mode characteristic impedance in the 
limiting case as the separation between the slots becomes 
very small. One expects that the even mode characteristic 
where Z, is 


impedance Z will be close to one half of Z 


0e (04 0 
the characteristic impedance of a single slot whose width 


is twice the width of the slot in the coupled structure, 
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Figure 2.2/2b Even and Odd Mode Characteristic Impedances 
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Figure 2.2/3 Dispersion Characteristics and Characteristic 


Impedance for Coplanar Waveguide with W/D = 1 
and ЕЕ 1.0 





The dispersion characteristic for both structures, however, 
should be the same. This comparison for the impedances is 
indicated in Figure 2.2/2Ь апа again shows an excellent 
agreement. 

It is interesting to observe the coupling and decoupling 
between the waves in the two slots as the frequency varies 
for larger valves of the separation constant S/D. With 
increasing frequency the waves become more closely bound 
to the slot which means that there is less interaction 


between the two waves. In this case Zoe and 2 converge 


00 
to Lor the characteristic impedance of a single slot with 
no coupling effect. 

Another interesting phenomenon is the fact that for a 
fixed КЕ the wavelength ratio À'/À in the even mode 
first increases and then decreases as S/D is increesed 
from very small to larger values as shown in Figure 2.2/2a. 
An explanation may be given as follows. For small enough 
values of S/D the metal strip between the slots has little 
effect on the propagating wave and the wave propagates as 
if it were in a slot of width 2 . W/D+S/D. Increasing the 
Separation between the slots effectively increases the slot 
width so that the metal strip still has little effect and 
the wavelength ratio increases. As the separation continues 
to increase, the two waves start to decouple and behave more 
as two waves on two slot lines which will finally be totally 
decoupled. Each wave then propagates on a slot line with 


width W/D, hence A'/A decreases. 
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However, again an error due to the first order approxima- 
tion can be observed. Note that for larger values of S/D 
the wavelength ratios and the characteristic impedances for 
the two natural modes should merge evenly together towards 
the respective line parameters for a single slot line. It 
is observed, however, that for still distinct wavelength 
ratios for the two modes the characteristic impedances for 
these modes merge not only together but overcross each other 
by a small amount. This fact is not shown in Figures 
2.2/2a-b for clarity purposes but is presented in Figure 2.2/4 
where for one arbitrarily chosen value of normalized frequency 
and slot width the two impedances for the two modes are 
displayed as a function of separation. Obviously the error 
introduced is not very large in this example but tends to 
increase somewhat for higher frequencies. Since the numerical 
methods used in this analysis were carefully checked it is 
believed that this error is caused by neglecting the 
z-directed component of the electric field in the power 
осма топ. As noted in Section 2.1.1 this component 
increases with frequency and the approximation caused a 
slight error in the power calculation cn a single slot line 
EwSection 2.1.2. 

However the partial comparison of the results with other 
data given above and an experimental verification described 
ШІ/Бессіоп 4 confirms the accuracy of the first order 


approximation for the frequency ranges considered here. 
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) COPLANAR STRIP LINE 

The configuration of a coplanar strip line is shown in 
Figure 1.1/4 and consists of two metal strip conductors 
placed on one side of the dielectric substrate. The spacing 
between these two strips is assumed to be constant in the 
z direction. The permittivity of tne substrate confines 
the fields to regions around the conductors to minimize 
radiation losses and to make this line Suitable for 
integrated circuit applications. A hybrid mode analysis in 
the spectral domain is again used to obtain the two important 
line parameters. However the key for a solution by this method 
is that the investigator has to have some idea of how the 
fields on this structure behave. In the study of the slot 
line an approximation ethe electr mer ricidmecr Gages НЕВЕ OE 
was required and was not difficult to obtain by phvsical 
reasoning. 

However for the coplanar Strip line an approximation of 
the current density across each strip is more feasible. 
Hence a set of coupled equations similar to (2-16) is 
desirable with the roles of current density and electric 
field interchanged. Note that this inversion was already 
obtained in the development of (2-16) and is the equation 


(A-22) in Appendix A give here as 


Гм. (а, в} Mo (a, B) J (a) E (a) 
T 2 


| 13 (0.8) м, (а, В) J (a) E, (0 
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Mecording to the scalar product definition and Galerkin's 
method dem in Section 2.1.1 a suitable set of weighting 
functions are the transforms of the current density distribu- 
tions. The methods of moments applied to (2-30) and an 
expansion of the current densities in basis functions yields 
the second order approximations similar to (2-19) 


со со 


с / мубо, в) |2,9) [аа +8 / M(a,8)3,(0)3,*(0)d0 


= СО 


N 
= 


(AC 


c S Mz(a,8)3,(0)3,*(a)da + a / M,(a,8)|7, (a “ао 


|| 
© 


= = СО 


In the detailed analysis of different approximations in 
Section 2.1.1 it was found that a first order approximation 
using only E, (a) was sufficient. Since the necessary power 
calculation for the computation of the characteristic 
impedance will only use a first order approximation, this 
approach is used throughout the analysis of coplanar strips. 
Thus J (x) ~ 0 and (2-31) reduces to the problem of finding 


 сорасабсіоп constant B such that 
f M,(G,B) |J, (o) |^àa = 0 (2-32) 
Bay One has the task here to choeoseva proper basis 
function. The following functions are reasonable basis 


functions for the current density component j, (x) over the 


strip conäuctors and have been used extensively in the 
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analysis of microstrip transmission lines. According tc 
Escure 1.1/4 the left strip extends for x such that 

ШИ - М < х < -S/2 and the right strip for 

5/2 <х < 5/2 + W respectively. Then the three possible 


basis functions are 


- left 
(17 (х) =1+1 right - 
Е, х) right (2-33a) 


0 elsewhere 








-i left 
W S+W, 2 
(5) = (х + i) 
(x) = E FOE (2-34a) 
(И) 2 "pr S+W, 2 
2 d 2 
0 еізеулеге 
= 
DOE oe 
-1 - W left 
2 | 
ES = (2-35a) 
MENU 
ШЕ 9 2 CGAC 
2 
0 elsewnere 


A Schematic presentation of these basis functions is given 


МР сиге 2.3/1. 
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Figure 2.3/1 


Assumed Current Density Component in z-Direction 
МШЕ ШШ с Сорїапагт Strip Line 

a) equation (2-33а) 

b) equation (2-34а) 

C) equation 2~35a) 
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The Fourier transforms can now be easily found by using 
again the shifting property cf the Fourier transformation 


indicated in equations (2-28). The transforms are given as 


(1) T P S+W E р 
Ја (a) = j a sin Ta у! . sin (<) (2-33b) 
TEE mus SH ui - 
I (a) mee sinia zn Bo | 57) (2-34b) 
2 
Ј03) (о) = 3 3 віп[а З] . [-24. + ХӘМ) sin 
(aw) ( aW) 
(01) 2—8 aW 
+ 3 со5(--=-) ] (2-355) 
3 2 
(aW) 


respectively. (Bg is the zero order Bessel function of the 
ЕЕ Капа.) 

It is obvious now that the theory and the numerical 
analysis of the slot line dispersion characteristic needs 
only minor changes to adapt a solution scheme for a coplanar 
Strip line dispersion characteristic. 

me compute the characteristic impedance of this trans- 
mission line a reformulation of (2-24) is necessary since 
the voltage between the two strips is not explicitly known. 
Again this definition is somewhat arbitrary because of the 
non-TEM nature of the problem. 


Define 
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where 1, is the total current on one strip and can be found 
analytically through an integration of the current densities 


(2-33a) through (2-35a) with respect to the variable x, 


yielding 
E = (2-37) 
CN 
Io = (a Sr 
160) 5. E 
Io E (2-39) 


The computation of the time average power flow follows in 
general the method outlined in Section 2.1.2. Some 
modifications are necessary, however, and are shown in 
detail in Appendix В. 

The results for the dispersion characteristic and 
characteristic impedance with these three different 
E mations are shown in Figure 2.3/2a-b. 

These line parameters are basically new. The wavelength 
ratio A'/\ was measured by Luna [Ref. 10] and good agreement 
is found between the present theory and the experimental 
data. Preliminary measurements on the characteristic 
impedance with a time domain reflectometer indicated good 
agreement with the thecry using the current density distribu- 
tion (2-34a). This work should be extended in the future. 

Reasonable agreement for the impedances may be found by 
comparing the present values with the quasi-static results 


by Wen [Ref. 4]. As one might expect, however, the present 
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0.46 4 
— _— еа  (2-33a) 
eq. (2-340) 
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method yields somewhat larger values due to the finite 
thickness of the dielectric substrate. 
2.4 PERTURBATIONAL ANALYSIS OF A SINGLE SLOT LINE 

ON A FERRITE SUBSTRATE 

The hybrid mode analysis presented so far has been applied 
to coplanar transmission lines on a dielectric substrate. 
Since these lines are new and only one complete frequency 
dependent analysis has been proposed it is not surprising 
that almost no analysis has apveared yet for coplanar 
transmission lines on ferrite substrates. 

The reasons why these lines offer interesting properties 
are Obvious and will be shortly given here. Ferrites are 
magetic oxides with high resistivity, a scalar dielectric 
constant anda tensor permeability at microwave frequencies. 
Electromagnetic fields within these ferrites will interact 
with the precession of electron spins, where this precession 
is related to the magetization vector which also precesses 
about the constant biaSing magnetic field. The magnetic bias 
field determines the orientation and the frequency of the 
precession and thus can be adjusted to any desired microwave 
frequency. Depending on the relative polarization and 
orientation of the RF magnetic field and the magnetization 
precession to each other, stronger or weaker interaction 
is obtained. Since the amplitude of the precession decreases 
due to damping the amount of energy absorption can be quite 


large for resonance. A Simple example of the use of this 
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phenomenon in ferrites is the isolator. If the RF magnetic 
field is circularly polarized and has the same orientation 
as the electron spin, strong interaction is obtained and a 
maximum of energy is absorbed. This attenuation would be a 
maximum for a wave travelling in one direction — where the 
orientation of RF field and electron spin is aligned — and 
smaller for a wave travelling in the opposite direction. 
Thus, the amount of absorption is dependent on the propagation 
direction which provides non reciprocal attenuation. Other 
non-reciprocal devices are switches, phase shifters, filters, 
and circulators [Ref. 19, 20]. 

It is then seen that coplanar transmission such as single 
Or coupled slot lines — or CPW — are bound to be used in 
ferrite applications since these lines posess a region of 
elanptically or circularly polarized magnetic RF field, much 
more pronounced than the quasi TEM microstrip structure. A 
Schematic representation of the fields on a slot line are 
shown in Figure 2.4/1. For these reasons and the fact that 
coplanar transmission lines offer advantages over the micro- 
Strip line as stated in Section 1, it can be foreseen that 
nonreciprocal devices for integrated circuits will be demanded. 
The need for an analysis of these structures becomes apparent. 
Furthermore, since the coplanar transmission lines as well 
as the properties of the ferrites are highly frequency depen- 
dent any analysis of these structures must be also frequency 


dependent. Іп this section an attempt is made to analyze the 
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slot line on a ferrite substrate by a perturbational method 
meme tne propagation cOnStant in the spectral domain. This 
method can then be easily extended for the analysis of 
coupled slots and coplanar Strip lines by the means described 
in Sections 2.1-2.3. 
СОО Perta on Expression forthe Propagation Constant 

Perturbational analvses are useful for the calcula- 
tion of changes in some quantity, if the problem itself is 
slightly changed or perturbed. Usually a solution for the 
unperturbed problem is known and the solution of the perturbed 
problem is found by assuming that these disturbances are 
small. Examples for the use of the perturbation theory are 
given in References 19 and 21. 

In the present application the unperturbed problem 
is the computation of the wavelength ratio A'/A — or the 
propagation constant — of the slot line on a dielectric sub- 
Strate. The solutions to this problem are known and given 
in Section 2.1.1. For the perturbed problem the dielectric 
substrate is replaced by a ferrite substrate with the same 
dielectric constant. It is assumed that this change in the 
Structure is small enough to allov a perturbational analysis 
of the change in the propagation constant. 

The perturbational expression for the propagation 
constant due to small changes in the type o£ material within 
the guidance structure is given by Helszajn in Reference 19 


as 


al 





-> 


іш SS a a a + (ро Aaa e Jas 
ү'+ү* = ы ————— е е (2-40) 


IJ (Ex х Н+ Е х H*).a, ds 
S 


where the following nomenclature is used 


S total cross section of waveguide 


AS cross section of waveguide which is perturbed by 
a change in type of material 


E and H are the electric and magnetic field of the unperturbed 
structure 


[AX ] change in magnetic susceptibility tensor 
[Axe] change in electric susceptibility tensor 
ү = 044438! propagation constant of perturbed problem 


ив X ) propagation constant of unperturbed problem. 


Although (2-40) applies in Reference 19 strictly to 
a closed boundary waveguide it can be shown that this 
expression also holds for an open boundary structure such 
as the slot line. This derivation is shown in Appendix C. 
A particular appiication of (2-40) for the slot line is now 
possible and is derived in the next section. 


28.2 Computation of Perturbation Expression 


me ee Sr em ee ee ee ee 


From the field configurations on the siot Jine it is 
seen that the magnetic field is elliptically polarized within 
the y-z plane. The strongest interaction between the RF 
magnetic field and the magnetization of the ferrite can then 


EN pected for a magnetic biasing field in tho x-direction. 





From the discussion of the susceptibility tensor in an 
ellipsoidal medium in Reference 19 it is seen that the 
EnSsent configuration is similar to that of a thin disk. 
Since the biasing magnetic field is along the x-axis no 
demagnetization for this field has to be considered. 

However, a demagnetization for the y-component of the RF 
magnetic field is to be noted. This was the deciding factor 
in using the external rather than internal susceptibility 
tensor, where the demagnetization enters into the formulation 
of the tensor rather than the fields and was found to be much 
easier to use in the following derivation. According to the 
small signal approximation of the equation of motion in 


Reference 19 two radian frequencies can be defined where 


WO = wr, ~ Nwo t КМ ш 
У 0 y 
OW = -N GU. + N_ uy 
Z 0 scm ZR 
Мс 
with wg 7 7Y i ; Oe Sey Yo 


у = -1.4 - g - 10° [Hz/0e] 
g Landé factor 
H 


magnetic biasing field in x direction 


XO 
Mo saturation magnetization 
N demagnetization factor (in this application DES N, = 0, 
ШЕ - 1) ; 
Y 





Furthermore, define the ellipticity of the normal modes of 


о) 
a e ze 
the uniform precession as eG =\ m and note that 
ct • ” \ M • * 
Dno ac Is#ätljeskittel resonancewrelatjione A@eeording 


to the Landau-Lifsnitz phenomenological damping formulation 


Y LH 
2 


DEO AH is the linewidth of the ferrite material. The 


[Refs. 19,22] the damping can be defined as au = - 





external susceptibiiity tensor is then 

















[x] = О Хе Xyz (2-41) 
E X2y X22 | 
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Equation (2-40) can now be evaluated. Note that 


the denominator can be expressed as EAM (see also 
fection 2.1.2) and that [Ах] = 0 due to the assumption 
that the dielectric constant remains unchanged, [Ax] 


is obviously identical to equation (2-41). To evaluate the 
remainder of the numerator ncte that the product 

[Ax] - H- H* has to be integrated іп the crosssectionai area 
LE ich extends such that -е <х <е апа О <у<р 

The RF magnetic field is given in terms of the scalar 
potentials in equations (2-2), (2-5) and (2-6). With these 


substitutions (2-40) becornes 





ju Uo =D д. оф 2 
1 AA Эз ДЕ ЫЙ... - Ма 
Ы” * Pavg s ) y P y E OI 
2 | aw ob, 
+ Е (+960 Jy T JWE- o 
6 2 P 
TT ! 1 
Ш X225 c2 527 I 
(2-42) 


The upper and lower signs denote waves traveling in the 
negative and positive z-direction respectively. 

This double integral is transformed via Parseval's 
relation [Ref. 15] into the spectral domain, the integration 
with respect to y computed analytically and the remaining 
Single integration is performed numerically on a digitali 


computer. Here again use is made of the íact that the 
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transforms of the scalar potentials are known after one has 
шта for the propagation constant B according to Section 
2.1.1. The time average power flow S is computed as 
explained in Section 2.1.2. Some details of the derivation 
leading to this single integral are shown in Appendix D, the 
method is the same as outlined in Sections 2.1 through 2.3. 

This perturbational analysis was carried out for an 
arbitrarily chosen set of parameters, such that resonance 
occurred at 4 GHz. Figure 2.4/2 shows the result of this 
calculation where B / B represents the ferrite phase constant 
normalized by the pure dielectric phase constant. The 
magnetic biasing field of Н — 683.6 Oe was directed into 
the positive x direction. The material constants are 
47M. = 2300 Gauss, AH = 100 Oe and E 2E 

This result in Figure 2.4/2 can be only qualitatively 
checked with results for similar structures for which a theory 
has been verified by experiment. Similar to the propagation 
constant for a ferrite filled parallel plane waveguide with 
transverse magnetization [Ref. 2] it is seen that there are 
two regions of propagation, below and above the resonance 
Frequency, separated by a region of high attenuation. It is 
noted that this attenuation is dependent on the direction of 
propagation. Another qualitative comparison can be made here 
with the resonance isolator described in Reference 19. Both 


Waves traveling in either direction are attenuated at 


4 


resonance, although one not as much as the other.  Accoràinc 
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to Helszajn [Ref. 19] optimum forward attenuaticn dces not 
correspond exactly to optimum reverse attenuation for 

E A En Since the normal modes are in general elliptically 
polarized. Due to the demagnetization described above this 
is the case here. Furthermore it must be noted that for the 
resonance isolator not to attenuate one of the propagating 
modes, the ferrite slab should be placed close to the er 
at some optimum position. But here the physical nature of 
the slot line determines an offcenter position of the ferrite 


if the conductor plane in y = D is thought of as a line of 


symmetry. 


2.5 SUMMARY 

A new method for a frequency dependent analysis c? open 
boundary coplanar transmission lines has been formulated. 
The two important line parameters, phase velocity and 
characteristic impedance, are derived from a complete hybrid 
mode analysis through scalar potential functions which 
satisfy Helmholtz's equation. This P.D.E. is solved by a 
Fourier transform technique and a matching of boundary 
conditions in terms of the electric field and current 
density components. The application of Galerkin's method in 
the spectral domain and an approximation of the electric 
field or current density by basis functions permits a solution 
for the phase velocity to any desired degree of accuracy. 
Through the application of Parseval's theorem to the integral 
of the complex Poynting vector in the space domain the 


evaluation of time average power flow on the waveguide 
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Structure in the spectral domain is feasible and allows the 
computation of the characteristic impedance of the line. 
The often cumbersome inversion of a Fourier transform is 
thus avoided. 

The only other existing analysis of a single slot line, 
Cohn [Ref. 3], has been numerically verified by a surprisingly 
close agreement between these two theories which are radically 
different. The present theory also provides new and more 
detailed information about the power density on the slot line 
Epucture. 

The formulation of the problem is so general that a 
simple change in the assumed field distribution vielded a 
complete analysis of two coupled slot lines. The data 
obtained is new and will be of major importance for the 
microwave industry in the design of slot line filters and 
couplers. Wherever a comparison with related theories was 
possible excellent agreement was found. Another modification 
of the existing theory provided the complete analysis of the 
coplanar strip transmission line. To some extent a comparison 
with experimental data was possible and snowed a very good 
Bereement. 

A perturbational spectral domain method for the analysis 
of the propagation constant of a coplanar transmission line 
on a ferrite substrate was developed and applied to the 
Specific example of a slot line. An experimental verification 


was not performed. 
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The spectral domain technique is relatively straight- 
forward in concept, but extensive algebraic manipulation is 
required to achieve computational efficiency. However, once 
the boundary value problems cf Helmholtz's equation is 
solved in terms of electric fields and current densities at 
the air-metal-dielectric interface, any coplanar waveguide 


structure can be easily analyzed. 
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2 m C T TOSED BOUNDARY STRUCTURES 


Lumped circuit elements are no longer of any practical 
use at microwave frequencies. Losses due to the increased 
resistance by the skin effect or due to radiation effectively 
terminates the use of Such elements. A common alternative is 
the construction of circuit elements by sections of trans- 
mission lines of which the quarter wave transformer is one 
of the most important examples. USing Standard transmission 
lines Such as coaxial cables or rectangular waveguide these 
lines are mostly short circuited at one end to keep radiation 
losses to a minimum. The short circuit can be assumed in 
most cases. to be ideal. This may be true and practical for 
closed boundary structures but is not easy to implement with 
an open boundary structure such as the slot-line. If a slot- 
line is simply terminated within the conducting plane by 
ending the slot it is found that a considerable amount of 
energy is stored behind this short circuit due to surface 
currents flowing around this end. This indicates that 
magnetic energy will be stored in this region so that the 
short circuit is not ideal but rather appears as an inductive 
reactance at a reference plane coincident with the end of the 
slot. This effect was first observed by Mariani and Agrios 
[Ref. 23] during their experimental investigation of slot 


line filters and later experimentally investigated by Knorr 
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and Saenz [Ref. 24] who also provided data curves of this 
so called "end-effect" for design purposes. No theory has 
appeared yet which describes this effect. 

Another example of how lumped circuits are replaced by 
alternative circuits at microwave frequencies is the 
resonating circuit.  Resonating structures are important in 
the design of filters and oscillators and have been constructed 
as cavities for rectangular or circular waveguides. For the 
slot-line to be an alternative transmission line to micro- 
strip line the use of resonating slots becomes evident for 
previously mentioned applications. Quarter wave coupled 
resonating slots and end coupled resonating slots have been 
presented on an experimental basis by Mariani and Aarios 
[Ref. DP put again no theory exists to provide means for 
the design. The problem here is how to determine the 
dimensions of a resonating slot such that resonance at some 
desired frequency occurs which is often the lowest possible 
frequency. Using the wavelength ratio of the slot line 
computed in Section 2 for this half wavelength resonator 
will not yield this information because of the above 
described end-effect. The wavelength ratio was derived for 
a transmission line, which extends infinitely in the 
z-direction (see Figure 1.1/2). 

An analysis of these kinds of problems is now proposed 
in this section. The hybrid mode analysis in the Fourier 


domain is used again as in Section 2 to provide a complete 
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frequency dependent solution. To minimize the losses due 
to radiation the resonating Slot section iS placed inside a 
shielding structure; this configuration is schematically 
shown in Figure 3.1/1. 

As a byproduct of this analysis the above cited end 
effect is solved theoretically. Some further modifications 
in the derived theory of the resonating slot structure 
unveils that line parameters Such as the dispersion 
characteristic and the characteristic impedance for shielded 
coplanar transmission lines can be obtained. AS an example 
results for a Single shielded slot line are provided. 

An important result of Section 2 is that once the boundary 
value problem of Helmholtz's equation is solved in terms of 
the electric fields and surface currents in the plane between 
the air-conductor-dielectric interface any similar structure 
can be easily analyzed by minor modifications. This result 
holds also for the closed boundary structures in this 
section and will be demonstrated with the example of end 
coupled resonating slot-lines. This structure can be 
described by a circuit model, where the elements of this 
model can be solved for. The result is a solution to an 
inductive coupling between two serial slot-lines. 

Throughout the following analyses in this section it is 
assumed that the two rectangular waveguides formed by the 
Shielding structures have dimensions such that rectangular 


waveguide modes are kept below cutoff. The cutoff condition 
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Fi | 
gure 3.1/1 Shielded Slot Line Resonator 
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for the waveguide formed by region 1 in Figure 398A is 
well known. Besides, if dimensions hy and ho are kept such 


nat only a TE mit mode (E field in y direction) could 


10 
be excited, the lower waveguide formed by regions 2 and 3 
will have the lower cutoff frequency due to the dielectic 
loading and has to be considered only. A cutoff condition 
was derived by Harrington [Ref. 21] for instance for a 
partially dielectric-filled rectangular waveguide. А 


solution is required for a transcendental equation given 


according to Reference 21 for the lowest possible mode as 


о К хз 
= tan (k 2P) ns cR tan(k, jh.) (3-1) 
3 
з и mE = D 
where k о = ш Еу, (37) (3-2а) 
С: о с 
EIE et iem 


The subscripts 2 and 3 refer to the regions 2 and 3 in 
Figure 3.1/1. This equation (3-1) is given here only for 
completeness and is solved numerically vithin the computer 


program described in Section 5.2. 


3.1 SINGLE RESONANT SLOT LINE 
The schematic of the resonating structure to be analyzed 
is shown in Figure 3.1/1 where it is assumed that the 


shielding structure extends for z such that -œ < z < œ, 





The conducting material on top of the dielectric and the 
surrounding shield is assumed to be perfect and the dielectric 
material is lossless and isotropic. 

Similar to the beginning in Section 2.1 a full mode 
analysis will require two linearly independent scalar 
potential functions. The z directed electric and magnetic 
field components in terms of these potentials may then be 


written according to References 13 or 21 as 


2 


О 2 
Ej (X1Y 12) = > + к. ) 94 (х,у,2) (3-3) 
2 
— 9 2 | 


i=1,2,3 defining the spatial region according to Figure 


21/1 


Ё c Um _ 
k = К = К = Ww 105) and Ko = € “Ko 


The transverse field components of the hybrid mode can then 
be derived through an application of Maxwell's curl equations 


to (3-3) and (3-4). 





Then 
HA | д, 
E  xcz my кик 
ИХ жи | 
Ее ______- = Bet ee e 
Eyi ду 92 “о Эх ыы 
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| 9%. дф, | 
Hyi 7 JU*i зу T Sx за aded 
| ДА ИХ 
Hyi = GE a + dy oz (3-6b) 


The scalar potential functions satisfy Helmholtz's equation 


in the three spatial regions [Ref. 21], thus 

= 0 en 
2 2 Е 

ТҮ ЕРКЕ 0 (3-8) 


where y? rare the three dimensional Laplacian operator 
wy, zZ coordinates., Two EE transforms are now used 
to transform these three-dimensional partial differential 
equations into ordinary differential equations which can be 
readily solved in terms of the boundary conditions. 

Define a two dimensional Fourier transform as 

o a ІШ ексе) 
G(a, ,y,£) = Л Л (хтио) е 2 dxdz (3-9) 
-9 -a 

Note that the transform with respect to x is a finite 
transform [Ref. 25], also widely known as a complex Fourier 
Series coefficient. 

Some physical reasoning must be used here to determine 
the discrete transform variable A, + For the lowest order 


mode xn this resonant slot structure the fields of the 





standing waves inside the slot in the z direction are similar 
to the z directed fields of a traveling wave shown in 

Figure 2.4/1. Thus, the z directed magnetic field component 
H, will be an even function with respect to x and from 
Section 2.1.1 and the second order approximation it is known 
that the z directed electric field component E, is an odd 
function with respect to x. A distinction in equation (3—9) 
in sine- and cosine-transforms aids in the following 


discussion. It is easily seen that if 


T 
Е(п) = J f(x) {ses (nx) dx ) 


sin 
-7 
then 
n r Z n 
J £"(x)cos (nx)dx = -n^F (n) + (-1) Е'(п) - f'(-1) 
-T 
and 
E 2 
f g"(x)sin(nx)dx = -n*G(n) + п[9(-т) - (-D'g(m)] 
-т 


So these transforms are expedient for problems in which 
either the function to be transformed or its derivative 

is prescribed at the endpoints of the interval. But note 
that for the above problem at the endpoints x = -a and x =a 
the tangential electric fields must be zero due to the 


conducting shield, hence E, i and E Or as Seen from 


yi' 
| дз , 
equations (3-3) and (3-5) А апа T vanish at these 
b ] ] * ~ • C . ~ • » 
oundaries since Ei 0i and H, i "E The transform 
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variable а, in (3-9) is thus determined to Бе 


With these considerations the two dimensional transform of 
equations (3-3) and (3-4) via (3-9) can be accomplished 
to obtain 


NES 
БОО У.Е) = (К, - ESO, la .Y,E) (3-10) 


H 73 (py Y, E) = (к; E xy 7. (а ,Y,8) (3-11) 


Similarly equations (3-5a) and (3-6a) are transformed into 


— 


oY. 
= æ == 7 __ 1 = 
Eug 711205 7 7616 5 (CE 
9%. 
H. = JWEs dy = A (3-13) 


бл 


The y directed components are of no concern in the development. 
Define 
2 2 20 2 | 
о, BE = ki ғуы (3-14) 


then the two Helmholtz's equations in (3-7) and (3-8) 


are transformed into 





ar. 2) $. (о E. | (3-15) 
қ 5 Hm 2 n'Y:$ 

Y 
22 ~ 7.7) v. (a ,y,E) - 0 (3-16) 
a 2 "ү + n'Y* 

у 
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It is observed, that the boundary conditions for the tangen- 


tial electric fields of y = D + h, and at y = -h, require 


1. 
E, = ( and Е, = 0, so that certain solutions to (3-15) and 
(3-16) are eliminated. These conditions are independent of 


x and z and may be applied in the transform domain as well. 


The solutions to (3-15) and (3-16) are then 


$, (o ,y, E) N sinh y, (D*h,-y) 


= е ] е ^ = 
Ф, (9 у, 5) = В (o. , 5) sinh Yoy + Cc (a, 18) cosh yoy (3 17) 


$3 (а ,y,E) = D’(a ,E) sinh y,(h,ty) 


— 


| h 
Y, (A, Y, 2) = А (a ,Е) cosh Y, (Dth,-y) 
i : h М 
Y, (ac ,y , &) = B (а ,Е) sinh Y Y + C (o, , E) cosh Y Y (3--18) 
h 
Y¿ (a, 1Y 18) = D (a +8) cosh Y, (hoty) 


It is to be noted at this point that n may assume the 


femme ZErO and -7% < E < о, d may therefore be less than 


2 


| 2 
zero for certain values of E. Furthermore, since к. > ку 


there will be three sets of solutions in equations (3-17) 


and (3-18) for n = 0, namely y, < 0, Y4^ « 0 where trigo- 


nometric functions replace all hyperbolic functions; Y, 220, 
2 | 
№ < Ó where trigonometric functions replace the hyperbolic 


functions in spatial region 2 and y, ко ү, > 0 то илеп 
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the solutions are as given above. It can furthermore be 

seen that these three types of solutions will exist only for 
п = 0. Assum п = 1 and E assumes its lowest value of E= 0. 
For s « 0 then Z < к. or “i < a. However, to avoid rectan- 
gular waveguide modes it is assumed that = See markus 

ram mO for n sl, £2, ... at alll frequencies of interest. 
According to the method of solution outlined in Section 2 the 
next step is the application of the continuity or boundary 
conditions at the interfaces between the three spatial regions. 
These conditions are given in equations (2-14) and (2-15) 

and apply for this closed boundary structure as well with a 
minor modification. The last four equations of (2-15) must 
now include a general dependence on the variable z for the 
electric field and current density components. A z dependence 
of e'f is sufficient for a ада ја о поли ine only. For 
instance, the current density component in the x direction 
would now be of the form jJ, 0062), the other components are 


changed Similarly in (2-15). Denote the transforms of these 


pouponents at y - D as 


Ela LE) E FiE, (x,2)) E, (a E) = F(E, (x,2)] 
Jy, (a, ' 5) = F5, (,2)7 J, (a, , E) = F(3,(x,2)) (3-19) 


If these components and the transforms of (2-14) and (2-15) 


are used with the substitution of the scalar potential 
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transforms through the solutions to Helmholtz's equations 
(3-17) and (3-18) it is possible to eliminate again the 
unknown coefficients AS through Dae Details of this exten- 
Sive algebraic manipulation are shown in Appendix E, the 


mesult of which is 


M4 (a, , £) M^ (a, , £) E, (0, 1 E) 2709,5) 
= (3-20) 
M, (a. , £) M¿ (a 18) E, (a, 18) Ј (о, 75) 
To apply here again the method of moments by Galerkin's 
method [Ref. 16] in the spectral domain the following 
definition of the scalar product is used [Ref. 15] 
сЕ(а '5);9(% 5) = ae Ты Ғ(а ,5)4%(а ,5) dé 
(3-21) 


For the application of Galerkin's method the suitable set of 
weighting functions is (E. (a, ,£) ,E, (a sd). With this set 
and through equations (3-21), (3-20) becomes 
( An ( 
О EJE, (a rEl Elo ,Ey» + <M la «EJE la +E)+E,la, rE)> 


(9-22) 


«M4 (a ,E)E, (a, ,E) ,E, (a, E)? i «M, (a, E) E, (a, 4E) , E, (n, , &)? 





The transforms of the field components are not known but 


can be expanded in terms of basis functions Such that 


Ela, = > а, Ё. ‚(а ,&) (52222) 
i-l 
Б (о 75) = Z b, ЕЁ ‚(о ,&) (3-23b) 
1=1 
Equation (3-22) can then be written explicitly as 
со со со E 
E Са - M) (og (E) Eu, (On EVE (0 15) 96 
со со со 3 
i s 2 o aT Mane Egi Onr E) Em OnE) dé = 0 
Nz= 2a (3-24) 
со со со У, 
E а : Ј M4 (a E) E i (a, ,E) E, (a, , E) dE 
ісі “ п--о -» 
+ т mE f Ma Gg E ESL 0S E) E, (91,5) dé = 0 


Equation (3-24) is still an exact formulation of the resonant 
problem to which a solution could be found in the following 
fashion. The coefficients М, to M, are besides other factors 


functions of the operating frequency, while the geometry of 


the structure and therefore the slot length enters through the 
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basis functions into the calculations. Given a desired fre- 
quency of operation the slot length L is changed until the 
determinant of the coefficient matrix of (3-24) is equal to 
zero. This still is a rather theoretical way of solving this 
problem since in (3-24) two infinite sets of basis functions 
are implied. Any solution will involve the numerical evalua- 
tions of the various coefficients, the integrations and other 
algebraic manipulations so that an approximate solution seems 
desirable. In Section 2 it was found for the open boundary 
structure that good agreement between this theory and others 
or experiments could be obtained if a good choice for a 

one term expansion of the electric field components was made. 
If the slot is relatively narrow the experience with the 
different approximations in Section 2.1.1 suggests the z 
directed electric field component at y = D may be neglected 
while the x directed field component at y = D may be approxi- 


mated such that 





А КО А (2 $ 
Е (х,р,2) = Е. (х,2) x е (х,2) = е] (х) еуі (2) (9225) 
where E W 
——— Ip s 
W, 2 2 
(ж) = х 
е (1) (х) = : (3-26a) 
xl / 
0 elsewhere 
1 TZ L 
к. ЕЕ Am |] < = 
е,1 (2) = (3-26b) 
0 elsewhere 
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The Fourier transform of (3-26a) was already given in 
equation (2-23a) but is presented here for completeness for 
ШИЕ finite transform. The transform of (3-26b) can be found 
Bnatytically. 


Ehe transforms of (3-26) are then 


(1) а n 
Е 1 (a) = T В (—5— (3-27а) 
EE 
T cos (== 
EN de) = (3-27b) 
х1 & 
2 2 
п’ - (EL) 


where Bo denotes the zero order Bessel function of the first 
kind. И 

It is бы possible to obtain a solution for the resonant 
slot length by varying L such that 


co 


г 7 мубо) [Е 01) (о) -Е(27 8) |2 ав = о (3-28) 


у == со = со 


This root finding can be accomplished now by a numerical 
solution on a digital computer, the details of which are 
given in Section 5. 

A representative graph of a result of this computation 
is shown in Figure 3.1/2 where the resonant length versus 
frequency for the lowest resonant mode is displayed. All 


quantities are scaled by the thickness D of the dielectric 
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Figure 3.1/2 Resonant Slot Length vs. Normalized Frequency 
for a Halfwave Slot Line Resonator with 
Y DA= 5.5, h,/D = h,/D = 8.0 and e, = 20 
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substrate which makes these results more versatile for any 
applications. These results are new and cannot be compared 
with any values resulting from other theories. However, 
values for the resonant length by this method are used to 
compute, theoretically, the end effect of a shorted slot 
line. Very good agreement with experiments is found which 
substantiates the validity of the theory and suffices as 


an experimental verification. 


Bee DISPERSION CHARACTERISTIC AND CHARACTERISTIC IMPEDANCE 

After the results for the resonating slots were obtained 
a thorough review of the involved formulas led to the 
following О dpa t which analyzes the line parameters 
for a shielded slot line. These parameters are of importance 
for mainly > reasons. 

The first reason is that for a theoretical analysis of 
the end effect of a single shorted slot-line the propagation 
Esnstant Of this line must be known. But so far only the 
open boundary slot line has been analyzed. An approximation 
could certainly be made in assuming that the boundary walls 
are far enough removed so that their effect on the slot line 
is negligible allowing the calculation of the end effect with 
the open boundary propagation constant. But for reasons of 
completeness line parameters for the shielded line should be 
EUolved to provide for studies on the effect of shields on 


the line parameters. 
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This effect now provides the second reason. If the 
slot-line is to be used in integrated microwave circuits 
physical dimensions will in general be small. Cohn noted 
in Reference 3 that the permittivity of the dielectric 
substrate must be sufficiently large to confine the wave to 
the slot and to minimize radiation losses. However, from 
the graph of the power density in the neighborhood of the 
slot shown in Figure 2.1/4b it is seen that even for = 16 
a substantial Cone of power 1s found for distances of some 
multiples of slot width away from the center. Hence, it can 
be anticipated that a shielding structure around the slot 
line seems desirable in some applications to comply with the 
Miniaturization of microwave circuits. Since no analysis 
of a Е ава slot-line has appeared yet in the literature 
the following study is another contribution towards a 
complete investigation of coplanar transmission lines. The 
method used follows exactly the one described in Sections 
КОШ and 2.1.2. 

A careful comparison of the equations used in Section 3.1 
for the resonant slot and those used in Section 2.1 for the 
Open boundary slot line reveals that only slight modifica- 
tions are necessary to adapt the theory in Section 3.1 for 
an analysis of shielded coplanar transmission lines, shown 
by the example of a single slot-line. 

If the general scalar potential function ф(х,у,2) 
in equation (3-3) is replaced by ф(х,ује Y? which accounts 


for a traveling wave in the negative z direction it is found 
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that equations (2-1) to (2-8) are equivalent to (3-3) to 
(3-8) if the partial differentiation with respect to z is 
accomplished, provided a similar replacement for v(x,y,2) 
of course. Now define the finite Fourier transform 

a о. © 


Р {ф(х,у)) = Ф(9 у) = S ф(х,у)е " ах (3-29) 
=d 


where o. is determined by the same considerations stated in 
Section 3.1. 
The Helmholtz's equations for the scalar potential 


functions are then 


Br 2 
I coa ) 2. (о, ,У) = 0 (320) 
у - 
al. - E ШЕ ) = 0 = 
5 ME i Се 
ду 
1c 
— 27 2 2 2 
where for the lossless case у = ЭВ and Ve IE em 


Hence equations (2-10) and (2-11) are obtained from (3-15) 
and (3-16) if Е is replaced by B. Similarly equations (3-10) 
through (3-13) hold for the transmission line problem if this 
replacement is carried out there. But then it is observed 
that an extensive algebraic manipulation towards a solution 
like that of Section 2.1 can be avoided since all equations 
Келатса in Section 3.1 up to equation (3-20) and all of 


Appendix E hold. 
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The final part of the formulation is then obvious. To 
solve for ß Galerkin's method is applied again, which 


requires the scalar product to be defined as 
== e * = 
<Е (а) ,9 (а) > 2 f(o)'g*(a.) . (3-32) 


As in Section 2.1, the electric field components at y = р 
are approximated by basis functions and used in connection 
with the above scalar product definition on equation (3-20) 
to obtain the set of simultaneous equations. The zero of 
the determinant of the associated coefficient matrix will 
provide the desired propagation constant ß. The details 
are by now familiar and given in Section 2.1.1. 

A first order approximation was used to obtain a solution 
where E, (x,D) zx 0 and E, (x,D) = E(x) was given by 
equation (2-22a). A representative result and comparison 
with an open boundary slot line is given at the end of this 
section. 

Inspired by the experience that only minor modifications 
were necessary to adapt equations from Section 3.1 for a 
Solution of the shielded slot line, the following investiga- 
tion is determined to evaluate the second line parameter, the 
characteristic impedance, whose definition is given in 
equation (2-24). The time average power flow on the shielded 


Slot line can be computed according to equation (2-26) as 
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h.+D a 
mer 1 IU 62 д, 2 900.2 9%, 2 
Pavg = 5 Ве = Е {онов СС) x (зу) ] zt oe [ (2) t 5) ] 
2 
Ты лл dy NE OS 23) 


ду 9х дх ду 


This formulation applies to the spatial domain and cannot be 
computed since the scalar potential functions are not known 
in this domain. Through the use of Parseval's relation 
(Ref. 15] it is possible, however, to transform this 
expression into the spectral domain where the Scalar potentials 
are known after the propagation constant 8 is solved for. 

Due to the finite transform in equation (3-29) the integra- 
tion with respect to x is transformed via Parseval's relation 


into an infinite summation, thus (3-33) becomes 





ds = Re к q (-a, ^Bue | 6^ - a, ^ Bug |? |” 
2 
= MES a ве |2514 
- Bl + 59 №9) 
+ ја К (0 E + = V)) dy. (3-34) 


According to the general procedure outlined in Section 
2.1.2 each coefficient of this series can be evaluated. The 


transforms of the scalar potential functions are given by 
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equations (3-17) and (3-18) where again the transform 

variable E is replaced by the propagation constant ß. The 
simple functional dependence on the variable y allows one 

to perform the integration in equation (3-34) analytically. 
Note that each coefficient of the series in (3-34) is actually 
the sum of three further coefficients which present the power 
carried in each of the three spatial regions. The details 
about these are shown in Appendix F. 

The final part in the evaluation of the characteristic 
impedance is shown in Section 2.1.2 and involves only the 
computation of the slot voltage to obtain the characteristic 
impedance according to equation (2-24). For programming 
purposes it should be noted that the computation for the 
shielded line is easier to perform Since a numerical integra- 
tion is avoided and replaced by a summation. The coefficients 
in the infinite series of equation (3-34) were found to decay 
rapidly so that a finite approximation yields good results. 

A representative set of curves from this calculation is 
displayed in Figure 3.2/1a-b where the wavelength ratio 1'/A 
and the characteristic impedance 20 of the shielded slot-line 
is compared with the line parameters of an equivalent open 
memmdary Slot line for two different slot widths. The cutoff 
frequency for the lowest order waveguide mode was computed 
according to equation (3-1) and is in this example in 
normalized form D/A = 0.0421. 


Several facts can be learned from this comparison. 
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It is observed that the difference between the parameters 
of open and shielded structures differ more for the lower 
frequency range than for the higher one, where the parameters 
of both lines merge together. This is no surprise if the 
following computation is performed. According to Cohn 
(Ref. 3] the ratio of the voltage V(r) along a semicircular 
path of constant radius г іп region l divided by the voltage 


V directly across the slot can be given as 





Mic) T (1) а 
vo = kor |H} (Кох) | (3-35) 
where E Cethe Rankel function of first kind and first 


order and 


е" Dus " Е 
ШІ - 727 5-5 cr) С) (3 36) 


Тһе Hankel function ни Ох asymptotically approaches 

е Х/ x for ~ > о so that an approximate result for the 
voltage ratio in equation (3-35) may be easily obtained. 

In the above stated example the waveguide walls are at 

r/D = 5,5 in the x direction. Then for W/D = 0.2 a voltace 
ШІП ОТ 0,34 іс computed for D/A = 0.02 and of 0.05 for 
ШІ! - 0.04. Similarly, the voltage ratio for W/D = 1.0 is 
0.39 and 0.07 for the normalized frequencies 0.02 and 0.04 
respectively at the position of the sidewalls. One may now 


conjecture that at the lower frequencies a considerable 
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amount of surface current flows on the sidewalls of the 
shields at |x/D|] = 5.5. This in turn means that due to 
these surface currents more power flow is distributed in 
the air regions 1 and 3 than is the case for the open 
boundary line. This conjecture is verified and can be seen 
in Figure 3.2/1 where the shielded slot line shows a larger 
wavelength ratio and a larger impedance for the lower 
frequency range due to this larger amount of power flow in 
air. The above calculation shows furthermore that the shields 
will effect a wider slot-line slightly more than a narrow one, 
@eeeto the closer confinement to the slot region of the latter. 
Again, this is verified by investigating Figure 3.2/1. 

As the frequency increases, the shields have less effect 
and the line parameters approach those for the open boundary 
line. However, the upper limit on the frequency of operation 


is the onset of a rectangular TE, 0 mode in the lower waveguide 


formed by regions 2 and 3. 


aoe oND EFFECT OF SHORTED SLOT LINE 

With the results of the preceding sections the end effect 
of a short circuited slot-line can now be computed analytically. 
A convenient description of this effect is based upon trans- 
mission line theory and can be performed in terms of a 
normalized reactance. This method was used by Knorr and 
Saenz [Ref. 24] in their experimental investigation of this 
effect and is accepted here to provide for a comparison 


between theory and experiment. 





Consider a slot-line which is shorted at one end by 
ending the slot in the conductor plane. The slot voltage 
would not only be zero at this short but also at multiples 
of A'/2 along the line if the short were perfect, thus 
forming a Standing wave pattern. Any section on this line 
between two voltage nulls could be considered then as being 
a half wave resonator described and analyzed in Section 3.1. 
If the half wavelength A'/2 of a wave on this line is compared 
with the resonating length L of the resonator for comparable 
physical dimensions at the same frequency, however, it is 
found that A'/2 > L. This could be expected since merely 
ending the slot in the plane of the conductor does not 
provide a perfect short circuit. Surface currents on the 
conductor beyond the extension of the slot-line increase the 
electrical length of the line so that the virtual short 
circuit appears at a distance A'/4-L/2 behind the physical 
short at the end of the line. 

The reactance caused by this effect is then computed as 
follows. The input impedance of a short circuited trans- 
mission line is 

Z. = j Za tan 6% (95:52) 
where 8 = 27/A' and ? = 21/4 - L/2. The normalized reactance 


can then be given as 


ГАЛЕ АА") т (3-38) 
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where the necessary data is provided by the analytic/numerical 
methods developed in Sections 3.1 and 3.2. Results of this 
calculation are shown in Figure 3.3/la-b and are compared 
with the experimental data from Reference 24. Although 

. actually the theoretical result for a shielded line is 
compared with the experiment from an open boundary line, 
excellent agreement can be seen. It should be noted, however, 
that the data presented was obtained by assuming that the 
z-directed electric field component could be neglected. It 
was found in Section 2.1.1 that this component becomes 
significant for normalized frequencies above 0.03. This fact 
may explain why there is a slight discrepancy between theory 
and experiment in Figure 3.3/lb for W/D = 0.541. A solution 
which accounts for the z-directed field component can be 
easily established since all necessary computations are 

given in Appendix E. This work was not performed here, 

Since the amount of computation time on the digital computer 
increases substantially. The good agreement between theory 
and experiment adequately verifies the method. 

The resulting data from this section are the basis for 
the theoretical design of microstrip to slot-line transitions 
which are based so far on the experimental data from Knorr 
and Saenz [Ref. 24]. Furthermore these results are used in 


Aec analysis of the following section. 
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БЛ END COUPEED SEOT-LINES 

An important results of Section 2 is that once the 
Helmholtz's equations in the various spatial regions are 
solved by the noted boundary matching technique in terms of 
the surface currents and electric field components in the 
plane y a D any configuration of conductors in this plane may 
be solved for the line or resonant parameters by approxi- 
mating either one of these field quantities. While in 
Section 2 this feature was demonstrated with the analysis 
of coupled slot lines and coplanar strip lines, the example 
presented here shows how still another slot-line configuration 
may be analyzed by a modification of the theory derived in 
Bection 3.1. 

The ES ding effect between transmission lines and 
resonators is used in the design of couplers and filters at 
microwave frequencies. This fact was mentioned befcre and 
is shown theoretically and experimentally in Section 4 for 
parallel coupling between two slot-lines. An analysis of 
parallel coupled slot line resonators for instance would 
basically follow the same procedure outlined in Section 2.2. 
Another configuration which has not yet been studied is 
investigated here: end coupled slot-lines. This end coupling 
effect between the two lines can be described easily in terms 
Of a mutual inductance since it was found in the previous 
Section that a simple short circuited slot line could be 


represented by an inductance at the end of the line. 
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The following analysis studies this coupling effect 
between two end coupled slot resonators in terms of odd 
and even mode resonant lengths. A Simple application of 
transmission line theory unveils the coupling effect in terms 
of network parameters such as end inductances and mutual 
inductance. A physical layout of two end coupled slot 
resonators is shown in Figure 3.4/la where the shielding 
structure is not displayed since the same configuration as 
ШШЕ иге 3.1/1 applies. A circuit model of the two slots 
with their mutual inductance can then be thought of as given 
in Figure 3.4/2. The analysis is then as follows. 

The impedance of an ideal transmission at any point on 


this line in terms of the load impedance Z the characteristic 


mn 
impedance 20 of the line and the electrical length fl is 


given as 


2, + 920 tan 82 


2 = 2 | 3-39 
0 20 + 92. tan BL 


where B is the propagation constant and £ is the length of 
the line from the load impedance to the point of interest. 
In this problem the load impedance Zr, is Z = Dr the short 
circuit reactance due to the end effect which is given as a 
normalized quantity in equation (3-38). Normalize (3-39) to 
obtain 

E ре 


7 = 


= Wc x, tan BZ (3-40) 
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а) 


b) 





Figure 3.4/1 End Coupled Siot Line Resonators 
a) Physical Layout and Dimensions 
p) Electric Field Component in x-Direction 
VS R Or Even and Odd Mode of Excitation 
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Since x tan (arctan хе) 


5 
tanlaretan x.) + tan 8% 
= = СЕЛІ 
x Па an x.)-tan BL 


or with the angie-sum relation of the tangent function 


р 
|| 


tan (arctan X, + 8%) (3-42) 


The network description of the circuit model in Figure 3.4/2 
is accomplished as follows. 

Using the two current loops in Figure 3,4/2 the corre- 
sponding loop eguations are 


— 


+ JoL I] * JoMI 


|| 
© 


41 


t 2 


(3-43) 


|| 
© 


ve 
21. + JwL I, + JwMi, 
men ter normalizing with respect to Zo: the characteristic 


impedance of the line 
(jx + Ix) Ty * jm I, = 0 


(3-44) 


tjml, о 
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Figure 3.4/2 Circuit Model for End Coupled Slot 
Resonators 


213 





Nontrivial solutions to this set of equations can be 


obtained for 


(x*x ) - m =0 
or 
А - = 
x +m e (3-45) 
Define LU (3-46a) 
ЕЕЕ (3-46) 


Then the inductances of each circuit loop and the mutual 


Enseetance of the coupling region are 


— 


ео Т 
ам 5 (3-47a) 
x = : 
and m = А (3-47) 


Equations (3-47) could be solved now if Хе апа хо аге known. 
These can be computed with the aid of equation (3-42) using 
the procedure outlined in Section 3.1 as follows. 

Let the even mode of operation be defined as 
E. (x, z) = Е,(х,-2) and similarly the odd mode as 
Е. (х,2) = Bat DE A plot of these functions is shown in 
Figure 3.4/lb. Given a set of physical parameters and a 


frequency of operation one computes the resonant lengths 


Lo and Lo for the even and odd mode of operation respectively 
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where the same method for obtaining these applies as outlined 
in Section 3.1 with equation (3-28) in general. However, due 
to the different geometry of this problem the approximating 
transforms of the electric field components in equation (3-27b) 


are now changed to 








(2) T cos = TAL 
E = 2 СО5(Е ) (3-46a) 
х1 12 Қ а)” 2 
EL 
(2) TACOS —— + 
Е 1 = ј 2 жо o sin (5 5 ) (3-46b) 
т" — (EL) 


for the even and odd mode respectively while equation (3-27a) 
remains unchanged. The shifting property of the Fourier 
transform was used as previously described in detail in 
Section 2.2. 

monce the two quantities L and Lj are solved for by a 
numerical root seeking procedure the even and odd mode 


reactances are then 


L 

a tan (27 Lt атса х.) (3-47а) 
Lo 

m tan (27 yr + arctan х.) (3-47b) 


The final computation for the reactances in the coupling 


region with equations (3-47) completes this analysis. 
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A representative set of curves resulting from this 
calculation is shown in Figure 3.4/3 where the following 
dimensions of the shielding structure were assumed: 

а/р = 5.5, h, /D = h,/D 20805 eU 20 and the normalized 
frequency of operation is D/A = 0.0159, which corresponds 

for the commercially available thickness of the dielectric 

D= 1/16" to a frequency of 3 GHz. As one expects, the 
normalized reactance in the coupling region is essentially 
meat computed for the end effect of acere re short circuited 
slot line, while the mutual inductance describing the coupling 
effect rapidly decreases as the separation T between the 

ends of the two lines increases. 

Although these results are not confirmed yet by experiment, 
ШЕ is Relieved that sufficient experimental verification 
exists in the comparisons of Section 3.3. The analysis of 


this section is nothing but a simple extension of the same 


method and the results should again be accurate. 


3.5 SUMMARY 

The theory developed in Section 2 has been extended for 
the analysis of shielded coplanar transmission lines and 
resonators. An exact formulation of the problem via a full 
or hybrid mode analysis is carried out in the Fourier 
transform domain where the transforms of dyadic Green's 
function components are derived. With the application of 
Galerkin's method in the transform domain and a finite 


approximation of the electric fields in the conductor plane 
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above the dielectric substrate a numerical solution to the 
slot line resonator problem is obtained. Although the 
primary intention was the analysis of a halfwave resonator an 
easy modification of the theory unveiled the line parameters 
of a shielded slot-line. Furthermore, with this data 
available the end effect of a single short circuited slot- 
line is theoretically determined. 

Once the Helmholtz's equations of the two scalar potential 
functions are solved for this inhomogeneous boundary value 
problem through a matching of boundary conditions in the 
various spatial domains in terms of the transforms of the 
current density or electric fields at the conductor plane 
above the dielectric any other configuration may be analyzed. 
The geometry of the problem encountered enters into the 
calculation only through the functional behavior of the 
assumed currents or fields. 

As an example of this result the coupling effect between 
two end coupled slot resonators was analyzed. Furthermore 
it is by now clear that any coplanar transmission line within 
the shielding structure may be investigated by simple 
modifications of the assumed distributions similar to the 
methods employed in Section 2 for coupled slot-lines and 
coplanar strip lines. 

The theoretical results of this section are new and 
provide important design information. An experimental 


verification is always desirable. Therefore, experimental 
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data for the end effect of a shorted open boundary slot line 
[Ref. 24] is used in a comparison with the theory developed 
here for a closed boundary structure. If the shielding 

walls are moved far enough from the slot it is intuitively 
clear that this comparison is valid. Excellent agreement is 
found for various parameters. Since the resonating length of 
the halfwave resonator as well as the wavelength ratio of the 
transmission line are involved in this computation, this 
experimental verification of the theory is quite thorough. 
Furthermore, shielded transmission line parameters tend to 
approach the values of the open boundary structure if the 
shields are further removed or similarly if the frequency of 
operation is increased. This also justifies the validity of 
the e ison. 

Since modification of this theory for the analysis of the 
end coupling effect between two slot lines is of a simple 
nature, further experimental verifications were not believed 
to be necessary. 

Unfortunately design curves for any of the configurations 
discussed in this section cannot be given for various reasons. 
Depending on a specific application the dimensions of the 
shielding structure can vary and hence will influence the 
line parameters. Line parameters are not only frequency 
dependent but depend on various considerations in specific 
applications, for instance in transitions to microstrip lines. 


So the parameters D/i, W/D, T/D and E, аге determined Бу 
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other constraints. Hence, the computer program proves to 
be more versatile than design curves and can be used in any 


specific application. 
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4A SLOT LINE DIRECTTONAL COUPLERIDESIGN 


a e Е-маил иде әс ы 


The design of passive components such as filters and 
couplers in the upper UHF and SHF frequency bands is based 
on electrical coupling between two or more transmission 
lines. (Some examples are bandpass filters, comb-line filters, 
interdigital filters and directional couplers.) The following 
analysis and experiments will be confined to the study of 
directional couplers, which are used as attenuators, power 
Splitters or hybrid junctions. However, most commonly 
directional couplers are used as a sampling device to measure 
separately the forward and backward waves on a transmission 
line, where a fraction of the energy traveling in either 
direction is coupled out such that the energy flow on the 
main line is not disturbed. Measurements of the reflection 
can be used in matching loads to a transmission line while 
the differences between forward and backward sampled power 
0а сасеѕ the net power transfer. The most often used 
Structure for the realization of a directional coupler is 
a four-part network formed by placing two transmission lines 
close enough to each other so that their fields interact 
[Ref. 18]. Numerous authors investigated these parallel 
line couplers with strip or microstrip lines. Although 
only recently the interest in the studies of these shifted 
from TEM analyses to frequency dependent solutions, one may 
say that the design and performance is well understood and 


known. 
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Since the introduction of the slot line as an alternative 
transmission line in 1969 by Cohn [Ref. 3] there has appeared 
Only one publication in the literäture about a slot line 
coupler. An experimental investigation by Mariani and Agrios 
[Ref. 23] points out that co-directional coupling can be 
achieved with two parallel slot lines. This fact is no 
longer surprising since Oliver in Reference 26 showed that 
contra-directional coupling is possible only for TEM lines, 
but co-directional coupling occurs if the phase constants 
for the two natural modes is different. That this is the 
case for the slot line is clearly seen now in Section 2.2. 

Since this thesis provides for the first time the 
analytical data necessary to describe the coupling effect 
between these two lines the objective of this section is not 
only a verification of the thoeretical results obtained 
in Section 2.2 by experiment but is also the first presenta- 
tion of a co-directional as well as a contra-directional 


Esupler using slot transmission lines. 


4.1 THEORETICAL METHODS FOR ANALYSIS 

A schematic diagram of the directional coupler to be 
analyzed is presented in Figure 4.1/1. Two of the terminals 
of this four-port network are in general connected to a main 
transmission line while the other two provide the sampled 
eueput. 

This four port may be described mathematically by either 


the scattering matrix or the impedance matrix where the 
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Figure 4.1/1 Schematic of Directional Coupler 
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elements of the matrices may be found from basic trans- 
mission line thecry. Both matrix formulations will require 
in general that the line parameters such as the character- 
istic impedances and the phase velocities for the two natural 
modes be known. Since in the design of directional couplers 
the ratios of power at the various ports are important to 
determine coupler parameters such as coupling coefficient and 
directivity, a more explicit formulation in terms of the 
ratios of the voltages at the various ports seems desirable. 
All four ports in general will be matched by the characteristic 
impedances of the connecting transmission lines, thus the 
load impedances for the four-port network are given by design 
Esucrderations and the ratios of power can be calculated from 
the voltage ratios. In the following sections two mathematical 
methods are derived, an impedance matrix method and a method 
which explicitly provides the voltage ratios at the various 
parts as well as the input impedance at port 1. The latter 
is used for computing the theoretical response of the coupler 
while the first method has been used aS a numerical check by 
an alternative method. 
4.1.1 Impedance Matrix Formulation 
According to Zysman and Johnson [Ref. 27] the impedance 


matrix description of the four-port network is 


И ТИЕ az) 





where [2] is a 4x4 matrix whose elements are given as 


231 = 222 = 253 7 244 И м сост) 
Z ей 7 = 2 = 2 = =j E (7 cot 0 - 7 cot 6 ) 
12 21 34 43 2 0е е 00 О 
2 = 2 = 2 = 2 = =) Er сәс 8» - 2 esc 5.) 
13 31 24 42 2 Ое е 00 О 
7 — 7 = 7 = 7 = =j B (Z csc 8 +7 csc 8 ) 
12 a 23 ma 2 0e e 0o O 


where 9. and ВЕ are the even and odd mode electrical lengths 
Sethe coupling region. 

Assume now that ports сс and 4 are terminated with 
some known impedance 22 and a known test voltage M is Supplied 
to port 1. Let the input admittance at port 1 be denoted by 
У. . Then a substitution of these values in (4-1) and a 


in 


reordering yields 





M лр E К q 7 
ШІ 1 2 25 22 іп 1 
2 7 7 
we 14 13 
44271 (te) E 5 0 
О О О 
= (4-2) 
7 7 7 
E. 11 12 
ds IER Gl Уз 0 
О О О 
2. 7 7 
BRE р 
Rom лен» ee | 67 | he 
О О О 
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ШЕ 1П а Shorthand notation 


ИЗИ [yj a) 


The unknown vector [x] can now be found by standard means 
for solving a set of simultaneous linear equations. It is 
seen that the matrices are in general complex. To avoid 
this complex arithmetic an isomorphism between this system 
of four complex unknowns and a system of eight real unknowns 
given by Zurmühl in Reference 28 was used. The imaginary 
quantities of (4-2) are then treated as real numbers. 
However, this method of describing the directional 
coupler is not very useful for design purposes, since not 
enough explicit information about the voltage ratios or the 
input impedance at port l is provided by analytical means. 
A more suitable formulation is given in the next section. 
Males Explicit Coupler Equations 
An explicit form of coupler equations was given by 
Jones and Bolljahn [Ref. 29], where the superposition prin- 
cipal was used. The input at port 1 is considered as the 
superposition of the two modes of excitation, odd and even. 
The analysis presented in Reference 29 is restricted in that 
it applies only for equal phase velocities of the two natural 
modes. By some algebraic manipulation it is possible to 
derive a more general form for unequal phase velocities. This 
calculation follows the general approach outlined in [Ref. 29] 


and details are not presented here. 


T26 





Assume that ports 2 to 4 are terminated with 





nt 


DE 2047206 (4-4) 


Then the ratio of the three voltages at these ports with 


respect to port l are given as 


lee 
8 cos Ве + J Б Sin ЫС | cos es HY R SIn 9 m 
Vi DE DO 
P. - - ер а 
Va DE DO 
1 DE DO 0 
where R= 2007 ° 08 


DE = 2cos I, u J (RIR) "Sin 9 


DO 2cos A в + 1/Е) sin 9 


The input impedance at port l is given by 


Уапи салон D j (B tan ош oC tan 9 о, 


2 іп E ^0 1 EE tan Б J = j(C tan s T B tan J с 
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where 





T aa ШЕШ ‚2 
22 
О 
в = 320e P 200 
20 
er E " 3200 
47 
О 


The necessary data for an evaluation of these equations is 
provided by the analysis in Section 2.2. The investigation 
Of Several configurations of directional couplers is now 


possible. 


meee THEORETICAL RESPONSE OF SEVERAL COUPLERS 

Some basic design considerations and the theoretical 
performance of several slot line couplers designed using 
data from Section 2.2 are discussed in this section. The 
analysis of the first coupler follows a similar investigation 
by Dalley [Ref. 30] for strip-line directional couplers. 
It is seen from equation (4-8) that a perfect impedance 
match at port 1 will be obtained if 05 = Zo: If any 
coupling between the lines is assumed then 2 де я 200 апа 
9 = 9 and an impedance match is Obtained only for 


tan 9. = tan S Or 
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A realistic example shows, however, that this design would 
be almost impractical even for n = 1 since a 10 dB coupler 
at 3 GHz should require a coupling length of approximately 
EE for E 16. Such a coupler shows the interesting 
property of power flowing mainly from port 1 to port 3 while 
a tenth of the input power is diverted to port 2; port 4 is 
essentially isolated at the center frequency. This 
theoretical nC is shown in Figure 4.2/1. It can also 
be seen that this coupler is much narrower in its bandwidth 
response than a conventional TEM line coupler. 

Another coupler Gesign could be based on the fact that 
the phase velocities of the two natural modes do not differ 
very inuch. An examination of equations (4-5) to (4-8) with 


~ 


the assumption 9. "р =з г тошо, = 0, unveils that one 
could indeed expect a similar performanc for this coupler 
as for a contra-directional TEM-line coupler, though not an 
Шоста се directivity. 

A digital computer simulation was performed and the 
characteristics of this simulated coupler are shown in 
Figure 4.2/2. This coupler has a wider bandwidth character- 
istic than the previous one and exhibits directivity of better 
than 6 dB. Due to the loose coupling the assumption 
20 = Ve T is still valid and the mismatch at port 1 
causes a VSWR of less than 1.08. 

Since the objective of this study is an experimental 


verification rather than an optimum coupler design no attempt 


to obtain a perfect match was made. 
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A third possible choice of design is one where 


B. T > = 21. This coupler will perform as a co-directional 


coupler with high directivity as can be seen if the 


following approximations are made. Assume loose coupling, 


then or AZ Zoo => with Zoe ~ 200° 


Е + а= 27 апа 9. is slightly larger than m while 9. 


is slightly less than т. Let 3 = sta and A = п-а 
where a is the difference between the electrical lengths for 
either even or odd mode and the electrical length of a 
single slot line. The physical length of the coupling region 
iS still computed from 8 Do s Ze 

Equation (4-6) can then be Simplified to 


V 


V 
> = a and equation (4-5) to E Е 0 The theoretical 
1 1 














performance of this coupler is shown in Figure 4.3/2. 


4.3 DESIGN OF TWO COUPLERS 

A comparison between the previously derived theory of 
a dispersive line coupler and its actual performance was 
made by designing two slot line couplers. A general layout 
for a coupler is shown in Figure 4.3/1. It was decided to 
use a microstrip to slot line transition to couple energy 
into and out of the coupler. The theory of this transition 
is well understood because Knorr in Reference 5 derived Я 
circuit model for this transition and developed a computer 
program by which an interative scheme vields m values of 


the various parameters required to achieve minimum VSWR. 
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Figure 4.3/1 Geometry of Experimental Slot Line Coupler 
with Microstrip-Slot Transitions 





This program was readily available and used in this study. 
The theory is aiven in detail in Reference 5 and therefore 
not reiterated here. 

The following design steps were taken before the actual 
physical construction could be made. 

a) It was decided to use standard miniature coaxial 
cable with an impedance of 50 Ohm to connect the coupler 
to the test equipment. This implied the use of a 50 Ohm 
microstrip line, thus the width of the microstrip was 
determined. 

b) The iterative procedure described in Reference 5 was 
used to determine the physical dimensions of the slot and 
microstrip line stubs and the width of the slot line at the 
desired operating frequency of 3 GHz for lowest VSWR. 

c) With the known dimension of the slot line the digital 
computer program described in Section 5.2.1 computed the line 
characteristics of two coupled slot lines for varying spacing 
(varying coupling) between the two lines. 

d) The data from step c) were used in another computer 
program which evaluated the performance of the slot line 
coupler but did not include the microstrip to slot line 
transitions. The coupling length could be varied. Both 
Mathematical formulations of Sections 4.1.1 and 4.1.2 were 
used in this design step. 

e) After the desired theoretical performance was obtained 


the lines were photo-etched onto the dielectric substrate. 
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Several iterations through Steps c) and d) are in general 
necessary until the requirements for the coupler are met. 

In the two experiments the coupled slots were etched 
onto one side of a l oz. copper surface on a dielectric 
substrate with thickness D - 0.125" and permittivity E 16, 
dimensions 4,5" x 6.5", while the four microstrip lines were 
etched onto the opposite side, to make all four ports 
accessible. OSM 248-4 adapters performed the transition from 
Microstrip to the miniature coaxial cable of the test 
equipment. A network analyzer was used for the recording of 
the coupler characteristic. 

The first coupler designed was the co-directional coupler 
shown in Figure 4.3/2. The discrete points in this graph 
indicate the experimentally measured coupler resronse. The 
slot width was 0.0625" and the separation between the parallel 
ЕНЕ 5 0,135". It was originally thought that the coupling 
length could be defined to be the length between the points 
A' and A" shown in Figure 4.3/1. This length was 0.375" and 
was found not to match the experimental with the theoretical 
results. It was learned, however, that good agreement could 
be obtained if the coupling length was taken between points 
вала В", in this case 0.75" which corresponds to 
m t 9% = 21. The waves on the two transmission lines do 
obviously not decouple abruptly in this particular design 
so that the coupling region is somewhat extended to an 
"effective" coupling length. The sidearms to ports 2 and 3, 
however, could not be bent further or the waves on these arms 


Nalda interfere with each other. 
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The second experiment was the design of a contra- 
directional 10 dB coupler. The slot width was again 
0.0625", but the separation between the slots was decreased 
22020375", The coupling length was taken between points B' 
and B" to be 0.371". The result of this experiment, shown 
by the discrete points in Figure 4.3/3 revealed some 
disagreement with the theoretical performance in Figure 
4.2/2. Although the experimental curves include the micro- 
strip to slot line transition, these transitions had a 
return loss of better than 10 dB over most of the frequency 
band and could not have caused this disagreement. Again it 
was found that the critical design parameter was the coupling 
length. Due to the closer coupling in this second coupler 
design, the coupling length had to be extended by 20% beyond 
the length B'- B" to match experiment and theory. Figure 
4.3/3 shows the comparison between the experimental and the 
adjusted theoretical response, where now better agreement is 
achieved. 

Although initially intended, design curves for slot 
line couplers cannot be presented since the slot line, 
single or coupled configuration, is a dispersive transmission 
line with frequency dependent characteristics. Thus any set 
of design parameters would be valid for a narrow frequency 


Band only. 





IN DB 


COUPLING 





2.0 3.0 4.0 
FREQUENCY IN GHZ 


Figure 4.3/3 Theoretical and Experimental Response 
of a Coupler with S/D = 0.3, W/D = 0.5 
and ne 16 


138 





4,4 SUMMARY 

An analytical evaluation of slot line directional 
couplers has been presented. Due to the different phase 
velocities of the even and odd modes, co- or contra- 
directional coupling can be achieved. Theoretical results 
were compared with experimental data by the design of two 
couplers. It was found that the coupling length is a critical 
parameter in this design. Coupled lines were assumed to 
extend infinitely in the computation of the line parameters 
in Section 2.2. However, in a coupler, abrupt decoupling 
See the two lines does not occur due to surface currents on 
the metal and this tends to increase the effective coupling 
length. This result is also new and will be important in 
future К бш. 

With an adjustment of the coupling length good agreement 
between the theoretical and experimental resuits further 


verifies the spectral domain technique used in Section 2 in 


general and the results of Section 2.2 in particular. 
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5. COMPUTER PROGRAMMING 


In Sections 2 and 3 the theoretical background of the 
analyses were provided together with some examples of the 
results obtained. By no means do these examples furnish 
enough information for any design purposes. Since this 
study presents the first extensive investigation of coplanar 
transmission lines, data aiding in the design are needed. 
Extensive design curves would not only be beyond the scope 
of this thesis but are even sometimes impossible to provide 
as stated before. For instance, for parallel and end coupled 
slot lines, not all applications can be foreseen at this 
Point. ПОС, a great effort was devoted to the developnent 
of two versatile computer programs which offer a variety of 
choices for the analysis of open and shielded boundary co- 
planar transmission lines. It is believed that these will 
be of more value to the microwave industry than any design 
curves could be. Before the two programs are described some 
general considerations which led to their development are 
necessary for the understanding of the numerical methods 


used. 


5.1 REMARKS ON NUMERICAL ANALYSIS 
The computation of the characteristic impedance and the 
perturbational analysis in Section 2 is based upon the solution 


to the dispersion characteristic problem for the transmission 


140 





line under consideration. In other words, the wave propaga- 
tion constant or wavelength ratio A'/A must Бе known before 
any other investigations can be started since only in this 
case are the scalar potential functions in the transform 
domain known. These are necessary for the evaluation of 
the power distribution and for the perturbational analysis 
for a ferrite substrate. Hence, in the analysis of slot lines 
the starting point is a solution for the propagation constant 
В to equation (2-19), while for copianar strip lines equation 
(2-32) applies. 

In general a 8 is sought such that 

со 2 
S Sta,ß) IT(a)| da = 0 (5-1) 
—оо 

where G(a,8) is a component of the transform of the dyadic 
Green's function N4 (a, 8) or M,(a,ß) with equations (2-19) 
and (2-32) respectively and T(a) 1S the Fourier transform 
of the electric field or current density assumed in the 
plane y = D for slot or coplanar Strip line respectively. 
A first order approximation is assumed here; if a second or 
higher order approximation is desired the zero of the deter- 
minant of the appropriate coefficient matrix must be found. 

It is obvious that an analytic solution in closed form 
for B cannot be obtained due to the algebraic complexity 
leading to the quantity G(a,8) as shown in Appendix A. Thus 
a numerical solution is mandatory and must be obtained by 


programming a digital computer. 





Several observations must be made to adapt the theoret- 
meal formulations in Section 2 to this numerical solution 
via a digital computer. 

It is shown in Appendix A that the transforms of the 
dyadic Green's function components М, - М, Or Ni - Ny are 
given in terms of coefficients F, through L, which are 
obtained by definitions during the application of the 
boundary matching technique. These coefficients are either 
real or imaginary so that in general a complex computation 
on the computer seems to be necessary. Furthermore, it is 


Een ain Section 2.1.1 that for all values of the transform 





арс jol < \x,2 - 8° the nature of the solutions in 
equations (2-12) and (2-13) changes from hyperbolic to 
trigonometric functions. For example, sinh YD changes 
to j sin Y5D . Since these functions enter into the 
computation of the quantities F, to N, it is found that these 
quantities are not only pure real or pure imaginary but 
change from real to imaginary or vice versa for the two 
cases Kee > 0 and ү,“ m on 

This fact furthermore seems to suggest a complex analysis 
on the computer. However, on the other hand, the disadvantaces 
of a complex arithmetic with digital computation are well 
known. Not only is the required amount of memory storage 
twice as large as for a real computation; a fundamental 
setback is that the computation times for instructions like 


mEiplications, divisions, and evaluations of trigonometric 


and hyperbolic functions with complex arguments are on the 
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average four to five times longer than in real arithmetic. 
An avoidance of digital complex computation is therefore 
highly desirable and can be obtained if some algebraic 
manipulations with the quantities F} to Na of Appendix A are 
performed prior to the computer programming. 

Before the details of this manipulation are shown a 
remark with regard to a scaling technique implied in the 
following derivation must be made. Dielectric substrates 
are commercially available for various thicknesses D so that 
in general the line parameters should be computed for each 
thickness and varying frequency. А more generalized method 
of analysis is simply to scale the various geometric 
quantities involved with respect to this thickness to obtain 
more ЖЕТЕ parameters for a varying normalized frequency 
D/A. The resulting design curves are then applicable for any 
geometry independent of the specific dielectric thickness. 

A manipulation of equations (A-13) to (A-24) including this 
scaling technique is given in Appendix G, where it is shown 
that the avoidance of complex arithmetic on the computer and 
thus long computation times can be traded against a somewhat 
lengthy but ые аи algebraic computation by hand 
before the ОО ои is attempted. 

With this formulation of the dyadic Green's function the 
further steps to obtain the wavelength ratio A'/A are as 
follows. For a first order approximation equation (5-1) has 


to be solved where now due to the scaling technique the 
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variable of integration is aD. The transform of the 

assumed electric field (or current density) for the geometry 
of transmission line to be analyzed is given in Section 2. 
Since in (5-1) the absolute value squared of these functions 
is to be taken, this part of the integrand will be in 
general a real and even function of aD. Preliminary numerical 
investigations showed that the components Mis Myr Ny and Na 
are also even functions of aD so that the whole integrand 

in (5-1) is even, a fact which allows further Savings in 
computation time since a numerical integration over a half 
interval between the integration limits is sufficient. Two 
representative examples of the real parts of the components 
M, and Ny and their dependence on aD are displayed in 
Figures 5.1/1 апа 5.1/2 respectively. Although these 
components according to the derivation in Appendix G are 
imaginary, this fact is immaterial if a first order 
approximation is used. 

For a second order approximation, however, the following 
observation must be made as explained here for instance with 
reference to equation (2-19). According to equations (2-21) 
or (2-23) E, (0) is an even and real function of a while 
E, (a) is odd and imaginary. N, and N, are found to be even 
functions while N, and N4 are ocd functions; all are 
imaginary. Equation (2-19) can be rewritten in a simplified 


form after the integration with respect to a is performed as 
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II о 
(5-2) 
A Piy = 
where the ШЕ are real quantities. A solution for the two 


coefficients a and b of the basis functions (cther than the 


trivial solution) does exist if п111022 = П]2°1>2] іп which 
T) 

case ја = Б M 
pan 


If b is assumed to be real this clearly indicates that 
a must be imaginary or vice versa, in other words, the electric 
field components in the x and z directions are 90? out of 
phase. Although this result is automatically implied through 
the notation used in the dispersion characteristic calculations 
it should receive special зо eatin in the power 
ealeulation. Since only a first order approximation for 
this computation iS-used in the present work, this phase 
shift was of no concern. 

A desired solution for equation (5-1) can now be found 
by assigning an arbitrary numerical value for the wavelength 
Шишо such that 1 < АМА! < с. and then iteratively 
changing this until the integration yields zero or more 
precisely a small enough value so that a prescribed accuracy 
for the wavelength ratio is realized. 

Several representative graphs of integrands for equation 
(5-1) are shown in Figures 5.1/3-6 where the coefíicients 


displayed in Figures 5.1/1-2 are used. Figure 5.1/3 shovs 
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the integrand for a single slot line, 5.1/4-5 two parallel 
coupled slot lines in the even and odd mode respectively, 

and 5.1/6 shows the integrand for a coplanar strip configura- 
tion. In all cases the zero order Bessel function was used 

as the basis for the transform of the assumed electric 

Mela (current density) distribution. 

The next task is the preparation of the appropriate 
equations for the numerical evaluation of the time average 
power flow. Here again complex arithmetic expressions can 
be avoided if manipulations similar to those in Appendix G 
are performed. General expressions for the time average power 
flow are given in Appendix B. Appendix H shows, how a real 
arithmetic evaluation of these expressions can be accomplished 
for the open boundary slot and coplanar strip lines. 

The normalizaticn of equations with respect to the 
thickness D of the dielectric (ferrite) substrate and the 
technique used to avoid complex arithmetic are given by пом, 
with sufficient explanation and detail so that no further 
expansion seems necessary. The perturbational expression 
for the ferrite analysis in Appendix D, for instance, need 
only be normalized. Appronriate remarks have already been 
made with regard to real and imaginary quantities involved. 
Hence equation (D-2) will require a real arithmetic only. 

Furthermore, although so far details are given for open 
boundary lines only, the manipulations described obviously 


apply for the shielded coplanar lines as well. The details 


To? 





of the preparation of these equations for programming are 
given in Appendix J. It is observed that for the dispersion 
characteristic computations the infinite integration in 
equation (5-1) is replaced by an infinite series due to the 
discrete nature of the transform variable O, > This remark 
holds also for the computation of time average power flow. 

Preliminary numerical investigations of the coefficients 
of these series indicated an even distribution with respect 
to the variable a, so that computation time can be saved by 
summing over a half interval only. 

This result even holds for computations concerning the 
length of a halfwave resonating slot, where equation (3-28) 
given in a simplified Form as 


со 


~ 2 
Br ER Му (9,5) IT (o4, £) | С 0 (5-3) 
has to be evaluated. Instead, due to the even behavior of 
the integrand with respect to 6 апа a. the following 
computation suffices: 


г у му (а, 5) [20а 5) [45 = 0 
nel. 0 


оо 


7 мусо, 5) |700,6) [2а + 2 
0 
(5-4) 


After this preparation of the equations for computer 
programming the only remaining task is to investigate the 
numerical integration with regard to its limits and the 


Procedure to be used. One of the limits of the integrals 





discussed so far extends to infinity sothat due to the 

finite presentation of numbers on a digital computer a finite 
approximation has to be used instead. From the representative 
examples shown in Figures 5.1/1-6 it is seen that the inte- 
grands will decay very rapidly in general. Although no 
analytic method can be applied to estimate the order of 
magnitude of a finite and adequate upper limit, a trial and 
error procedure by numerical means showed that for the open 
boundary lines a limit of aD = 55 and for the shielded lines 
ED = 9 is sufficiently large to be considered infinite for 
all practical purposes Since the integrands are at least down 


5 


by a factor of 10 ~ relative to their largest value close to 


the origin. Several different investigations showed that 
these binis guarantee at least a three digit accuracy for 
any practical set of line parameters. 

After establishing appropriate limits of integration 
the remaining problem of which numerical integration scheme 
to use must be considered. Many schemes are supplied through 
the standard literature on numerical analysis. In this study 
only three methods were investigated seriously. Two 
integration subroutines are supplied as part of 'IMSL' 
Library [Ref. 31] with names 'DRMBIU' and 'DCADRE' where 
both use a repeated interval-halving technique, also referred 
to as a Romberg integration.  'DCADRE' is an improved version 
of 'DRMBIU! since it provides for an adaptive method and 
thus may save computation time. This routine is quite new 


and became available only very recently. It is used for 


comparative purposes in this discussion. 
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An alternate integration scheme is the Gaussian- 
quadrature method. According to Reference 32 the Gauss- 
Legendre formulas are significantly more accurate than the 
equal-interval formulas such as the Romberg-method for a 
comparable number of sub-intervals between the integration 
limits with the condition that the magnitudes of higher order 
derivatives of the integrand decrease or at least do not 
increase. In other words, an equal-interval formula will 
in general need more divisions into sub-intervals to acquire 
the same accuracy which results in a greater computation 
time. Since the computer programs presented here are intended 
to be used very often for design purposes, the computation 
time is an important factor. Considerable effort was devoted 
therefore to finding the optimum numerical integraticn scheme 
where “optimum" means the shortest possible computation time 
for a prescribed accuracy of the result. 

Due to the algebraic complexity of ane integrands the 
behavior of its higher order derivatives cannot be studied 
on an analytical basis. It is observed, however, that these 
integrands will partly consist of the transforms of the 
dyadic Green's function components. The examples shown in 
megures 5.1/1 and 5.1/2 indicate a smooth and continuous 
behavior with at most one zero for oD in the interval [0, ). 
Bar zother part of the integrand will consist of the transform 
of the assumed field (current density) distribution. If, 


for instance, coupled slots in the odd mode are considered, 
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Еріс function will be in the form of 


S/D + W/D 


вїп (#^——————) ар) т ee 


5 ) . Let y = E(x) where y is 


the largest integer less than or equal to x. Then the first 
term above has E (22 (8/D +W/D)) zeros for Dain the interval 
(0,55]. The second term has approximately E (32 W/D) zeros 
Since for large arguments s, Bg (x) ~ cos(x). These zeros 
are not necessarily coincident with any of the first ones 
so that the total number of zeros of the integrand is the 
sum of both. 

For resonant slot length computations that part of the 
integrand involving the transform variable ED is proportional 


ED L/D 
(22.02) 


to cos (see equation (3-27b)), a function which 


I zeros for ED in the interval (0,9]. 


Since in general the integrands will be basically in the 


will have E( 


КООШО а product of a trigonometric апа а zero order Bessel 
function the following test was performed to find an efficient 
numerical integration scheme for the specific needs in this 
analysis. 


A test function to be integrated in the interval [0,50] 


was chosen to be f(x) = Sin (5x) *Bp (2x) which has 68 zeros 
Aithin this interval. f a three digit accuracy is desired 


a numerical Romberg integration with 'DRMBIU' needs 
approximately 6 seconds computation time while the cautious 
adaptive Romberg method in 'DCADRE' needs only 4 seconds. 

On the other hand, a Gaussian-quadrature 80 points integration 


requires only 0.5 seconds to guarantee the same result. 





The advantage of this integration scheme becomes obvious 
then. However, in contrast to the Romberg integrations, 
where the user can simply specify the accuracy of the result, 
this is not possible with a Gaussian-quadrature routine. 

Here the user must have advance knowledge of the general 
behavior of the integrand to select the number of quadrature 
ponts. Still, no estimate of the integration error can be 
obtained. 

But from the numerical example above it was conjectured, 
that in the present applications a quadrature integraticn of 
sufficient accuracy can be achieved, if the number of quadra- 
ture points is at least twice the number of zeros of the 
integrand within the limits of integration. 

This number depends on the parameters such as slot width, 
separation between slots or strips and for the resonator 
analysis on the frequency of operation and will usually vary 
over a wide range. ^A basic feature in the first program 
'COPLAN' is therefore a short computation which estimates 
the number of zeros of the integrand according to which the 
selection cf an appropriate quadrature routine is made. It 
was found by several test runs and comparisons that a sub- 
Stantial amount of computation time could be saved without 
any Sacrifice in accuracy. 

5.2 COMPUTER PROGRAM 'COPLAN' FOR OPEN BOUNDARY COPLANAR 
TRANSMISSION LINES 
A numerical analysis of open boundary coplanar trans- 


mission lines is provided by the digital computer program 
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'COPLAN' and itS accompanying subroutines ‘COMBO' and ‘'GAUSQ', 
all shown in Appendix M, written in the FORTRAN IV programming 
language. The Gaussian-quadrature integration routines are 
not shown. 

Although during the development of the theory in Section 
2 of this thesis separate programs were written for different 
line configurations, this program 'COPLAN' is a combination 
of all of them and offers a high degree in versatility.  Dis- 
persion characteristic and characteristic impedance may be 
Semputed for .Single or coupled slot lines, for coplanar 
waveguides or for coplanar strip lines. Except for the 
latter the perturbational analysis for phase shift and attenu- 
ation may be applied as well for ferrite substrates. 

The necessary input data and some further information 
about the use of this program are described in the preceeding 
commentary section of the main program. Input and output 
formats are chosen such that design curves may be easily 
generated if so desired. 

After the input data is read in and printed in an echo 
check, preliminary computations and initializations of the 
various computation loops are performed. The first computa- 
tion inside the inner loop is a check whether the operating 
frequency is above the cut-off-frequency of the lowest order 
surface wave mode; if so the job is terminated. The next 
step is to obtain an estimate of the number of zeros to be 


Eucountered in the interval of integration, according to 
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which a numerical value is assigned to an integer variable 
to be used in the selection o£ the appropriate Gaussian- 
quadrature integration routine. This selection and the 
actual integration are performed in the subfunction 'GAUSQ' 
whose use is self-evident. 

The program proceeds with a zero finding procedure to 
find the wavelength ratio A'/A, according to equations (2-19) 
and (2-31), where 1 < A/A' < Е. Integration limits are 
computed to determine for which values of aD the trigono- 
metric and hyperbolic functions in the solutions to Helm- 
holtz's equation apply (Section 2.1.1). For hyperbolic 
BUHCtions (K 5D < |ор| < 55) the appropriate numerical inte- 
. gration is performed via the subroutine 'GAUSQ' while for 


Шие |ор| < ko D it was found that a 24 point Gaussian-quadra- 


2 
ture routine is accurate enough for all practical parameters. 
The subroutine COMBO provides the appropriate integrads via 
several identifiers or flags. 

After the integration is carried out the result is used 
to select a new value for A'/A through a bisectional zero 
finding to improve the result and the procedure starts over 
again. The computation is terminated when the absolute value 
of the difference between the new and the old value for 
A/A' is less than O and thus a three digit accuracy 
after the decimal point is guaranteed. 

Preliminary computations follow for the determination 


of the characteristic impedance. Depending on the trans- 


mission line configuration, mode of operation and the spatial 





region the appropriate integrands are provided again through 
the subroutine 'COMBO' in accordance to the equations of 
Section 2. The numerical integration of the power density 

in the spectral domain follows the same considerations given 
Bbeve for the dispersion calculations. The proper definition 
of the characteristic impedance for the transmission line 

to be analyzed is then used to evaluate this line parameter. 
A three digit accuracy is quaranteed in this computation. 

The two line parameters for a dielectric substrate are 
thus determined and are printed into a tabulated output. If 
a perturbational analysis for a ferrite substrate is desired, 
however, the computer program proceeds with preliminary com- 
putations to evaluate the external susceptibility tensor 
components. These are then provided to the subroutine 
'COMBO' which computes the appropriate integrands for the 
perturbational analysis of the attenuation and the phase 
shift for waves travelling in either direction, After the 
numerical integration is performed some further manipulations 
are necessary to allow the attenuation to be printed in the 
unit [dB/cm] and the phase shift as a ratio of the ferrite 
phase constant over the dielectric phase constant. These 
results are printed together with the results for a dielectric 
Substrate. The procedure is then reiterated for the next 
set of parameters. 

The subroutine 'COMBO' provides the integrands for the 


dispersion characteristic and power density computations as 
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well as for the perturbation analysis. The entry variabie 
of this routine is the normalized transform variable aD 
while other coefficients are exchanged in common storage 
areas (COMMON statement) with the main program. Depending 
on the transmission line configuration the appropriate trans- 
form of the assumed electric field (current density) is 
computed by the use of a Bessel-function subroutine. The 
further evaluations of the respective integrands follovs the 
general theoretical outline of Section 2. Variable names 
are chosen such that these equations can be easily compared 
with the derivations of equations in the appropriate appen- 
dices so that no further explanations are necessary. 

Representative computer outputs are given in Tables I 
through VII. Waves travelling on a ferrite substrate show 
а попгесїргоса1 behavior. It is therefore pertinent to cite 
psexorrentation of the coordinate system and the biasing 
field. Following Оо ЕП сонное Section 2.4 the direction 
of the coordinates are given in Figure 1/2 and the biasing 
magnetic field is in the positive x-direction. A forward 
wave then travels into the positive z direction. 

Finally some statistics about this program seem to be 
apt. The total program including all subroutines requires 
182 K bytes storage space and 78 K bytes for the execution 
Step. To obtain a fast execution time the program was com- 
piled with the FORTRAN H compiler which applies an optimiza- 


tion code. The computation time for execution vary greatly 
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COUPLED SLOT LINES ON FERRITE SUBSTRATE 
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with the line configuration and the physical parameters so 
that only some typical examples are given here. To generate 
the Tables I through VII the respective computation times 
ШЕ 60, 65, 95, 45, 55, 300 and 145 seconds. 

For a comparison Of accuracy, storage and computation 
time the Gaussian-quadrature integration was replaced by the 
adaptive Romberg integration routine 'DCADRE'. The total 
program length was found to be 215 K bytes while the computa- 
tion time increased at least fivefold to guarantee the same 
numerical accuracy. The detailed and involved investigations 
of the numerical methods are thus justified since with the 
here presented program an efficient numerical application of 
the theory of Section 2 is presented and can be used Бу 
Future users. 

5.3 COMPUTER PROGRAM 'RESO' FOR RESONATING SLOT LINE AND 

TRANSMISSION LINE PARAMETERS FOR SHIELDED CONFIGURATIONS 

A computer program which provides a numerical analysis 
of shielded slot lines and halfwave slot resonators is shown 
in Appendix N with the main program called 'RESO' and its 
accompanying subroutines 'COMBC' and 'DET', all written in 
the FORTRAN IV programming language. The preceeding commen- 
tary section of the main program describes the use, the neces- 
sary input data and explains some of the variable names which 
are important for input/output operations. 

After the input data is read the program computes the 
lowest order mode cut-off frequency for a rectangular wave- 


guide mode in the guide formed by the shielding walls and 
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the metal above the dielectric substrate, thus including the 
latter. The determining equation is given by (3-1) and can 

be easily programmed, shown in the subroutine 'DET'. The 
result is printed and is subsequently compared with the 
Operating frequency where any further computation is terminated 
if the operating frequency can cause an onset of a rectangular 
waveguide mode. Initialization of computation loops for 
varying normalized frequency and slot width follow. 

The first major computation is to find the wavelength 
ratio A'/A. The necessary coefficients for the determining 
series are provided by the subroutine 'COMBC'. Similar to 
the procedure described in Section 5.2 a bisectional search 
to find the zero of the determining functional is used to 
improve subsequently the accuracy of the result through 
several iterations. The computation stops when an accuracy 
of three digits after the decimal point is guaranteed for 
the ratio A/A'. The program proceeds with the computation 
of the characteristic impedance. The coefficients for this 
computation (equation 3-34) are provided again by the 
Subroutine 'COMBC'. 

Finally the computation of a halfwave slot resonator is 
Started. Similar to the iterative procedure for the disper- 
sion characteristic, the halfwave resonating length is found 
by a bisectional zero finding procedure. The result from 
the dispersion characteristic is used to ensure that the 


Borrect zero is found since the slot resonator would not 
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only resonate at a halfwave length but also at multiples 


thereof. For a halfwave length resonator the normalized 


ATA 
2D/A 


improved until a one digit accuracy after the decimal point 


resonating length L/D < A first guess is subsequently 
for L/D is obtained. The necessary integrands for equation 
(3-28) are provided by the subroutine 'COMBC'. It was men- 
tioned in the general considerations to the numerical methods 
employed, that an estimate of the number of zeros of the 
integrands could improve the computational efficiency. By 
investigating various choices of input parameters it was 
found, however, that the following number of points of 
Gaussian quadrature integration schemes guarantee a three 


digit accuracy for the integration: 


y, «m, ү, < 0 2 Points 
2 2 ! 

Yi > 0, Yə < 0 24 pcints 
2 2 | 

were 0, e 0 56 points 


After finding the resonating length the program proceeds 

with the computation of the end effect of a simple short 
circuited slot line by evaluating the equations for the 
normalized short circuit reactance. All results are printed 
together with some of the input data. If no analysis of end 
coupled slot lines is desired the total procedure is reiterated 


for the next set of parameters of frequency and slot width. 
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The analysis for end coupled slot lines is performed in 
the same manner as the single slot resonator length is 
computed.  Normalized reactance and mutual reactance of the 
coupling region are derived from the resonating length in 
the odd and even excitation of two resonating slots, For a 
given operating frequency and slot width these parameters 
are computed for the desired number of different separations 
between slots and are printed. The program starts for 
another set of input data from the beginning. 

The subroutine 'COMBC' provides the coefficients and 
integrands for the various computations. The entry variable 
is BD for dispersion and impedance computation and is ED 
for the resonator analysis. All variable names have been 
chosen such that the coded statements can be readily compared 
with the respective equations of the appropriate appendices. 
The transform of the assumed EE T field distribution is 
chosen as follows. It was learned, that for the dispersion 


characteristic and resonator analysis a distribution 


3 
x 
1 + [475 |х| < м/2 


Жа E (5-5) 


0 elsewhere 
which is mathematically similar to (2-35a) gives almost iden- 


Ecol results as if a distribution according to (3-26a) is 


chosen; both analyses seem to be relatively insensitive to 
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the choice of basis functions while the characteristic 
impedance computation is more sensitive. Furthermore, 
observe that the largest computation time in this program is 
needed for the iterative procedures for dispersion and 
resonator analyses. The numerical evaluation of the trans- 
form of (3-26a) given in (3-27a) requires the use of 
another subroutine — a Bessel function routine — while the 
transform of (5-5) simply involves algebraic manipulations 
of trigonometric functions as can be seen in equation 2355) 
and thus is faster on the computer. This aspect is observed 
in the subroutine 'COMBC' which uses the Bessel function 
for the power density computation only. 

The total program length ans all subroutines is 
54 K bytes in storage. The program was compiled with the 
FORTRAN H compiler. A sample output is shown in Table VIII, 


where 23 minutes computation was needed. 
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6. SUGGESTIONS FOR FURTHER RESEARCH 


Although an extensive analysis of slot lines and coplanar 
strip lines has been presented, Several problems still remain 
which could be the subject of further investigation. The 
Intention at this point of the thesis is to simply list a 
few of these and to provide some guidelines for their 
possible analysis. 

The transforms of the dyadic Green's function components 
which were derived here are general and could be used in the 
analysis of different coplanar configurations. For instance, 
Mariani and Agrios presented in Reference 23 some preliminary 
experimental results for quarter wave coupled slot iine 
resonators in a shielding box. This configuration could be 
useful in the design of microwave filters. To obtain 
numerical results for this type of Mime all one has to do is 
to approximate the anticipated electric field distribution. 
This approximation has to be Fourier transformed and is then 
applied in equation (3-28) where it replaces the transform 
of the assumed electric field for a single halfwave slot 
line resonator. 

Coplanar transmission lines on a ferrite substrate have 
been analyzed in this thesis by a perturbational expression 
for the propagation constant only. These results still need 
an experimental verification. Furthermore, an exact analysis 


for these configurations seems highly desirable and is still 
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missing in the literature. It seems that this formulation 
can be easily obtained in the following manner. The 
transforms of the dyadic Green's function components should 
be derived for a ferrite substrate. The derivetion vould 
follow exactly the general concept shown in Appendix A where 
the permeability tensor characteristics would be implied in 
a general formulation from the start. Although this method 
would result in a higher degree of algebraic complexity it 
should be by now a straightforward procedure. This work 
would be more general and complete than any of the two 
analyses reported so far [Ref. ll, this thesis]. 

Circular apertures are used in RS engineering and 
their analyses are described by various authors in the 
ШЕ... But nothing has been reported so far about a 
erreular aperture On dielectric substrate апа цепсе 
specific properties are unknown. An attempt was made to 
apply the Fourier transform method derived in this thesis to 
an analysis of this configuration. The steps of this analysis 
are briefly outlined in Appendix K and indicate that various 
and difficult problems are encountered. It cannot be 
predicted whether this formulation is feasible at all. 

A final suggestion is the analysis of sandwich slot or 
coplanar strip lines. These line configurations will be 
К lar to the one described in this thesis where another 
electric slab is placed on top of the conductor. This 


configuration was first suggested by Mariani and Agrios 
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[Ref. 23) since it offers the advantages of shorter wavelength 
and greater confinement of the electromagnetic fields and 

thus can reduce radiation losses or perturbing effects from 
nearby obstacles. A possible analysis vould be started by 
developing the transforms of the dyadic Green's function 
components for this geometry where perhaps some simpiifica- 
tions could be made due to the symmetry with respect to the 


conductor plane. 





7. W@ONCLUSION 


A complete hybrid mode formulation has been presented 
for the analysis of coplanar transmission lines. The 
analytical methods are applied in tne Fourier transform 
domain where dyadic Green's function components have pecu 
derived for open boundary and shielded transmission lines. 
These components are very general and do not depend on the 
Specific geometry of any particular line configuration. 
Using these in connection with Galerkin's method, exact 
equations are obtained from which line parameters such as 
the length of a resonator and dispersion characteristics 
of transmission lines can be computed. Approximations by 
the use of a finite set of basis functions yield good 
numerical results. 

Using the results of the dispersion problem the theory 
Was extended to a formulation of the characteristic impedance 
in the Fourier transform domain. In still another extension 
a perturbational expression in the spectral domain was applied 
for the analysis of the propagation coefficient of coplanar 
transmission lines on a ferrite. substrate. Both of these 
applications show the flexibility cf the formulation in the 
Fourier domain since all necessary computations can be 
performed in this domain. Inversion of the transforms, a 
task which quite often is difficult to accomplish, is thus 


avoided. 
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Some design curves and two computer programs have been 
derived and are provided within this thesis. Various new 
information for the design of coplanar transmission lines 
is therefore available. 

Where applicable, results were compared with other 
existing theories and good agreement was ag all cases. 

Comparisons with experiments were also performed. In 
particular the numerical line parameters for coupled slot 
lines were used in the design of slot line directional 
couplers of which several different types were investigated 
theoretically. Due to the dispersive character of a slot 
line co- and contra-directional coupler design is feasible. 
Two directional couplers were built and their actual 
performances measured. Good agreement between theory and 
experiment further verified the validity of the analyses. 

In this thesis particular interest was devoted to the 
analysis of specific coplanar transmission lines such as 
slot lines and coplanar strip lines. However, the method 
employed is completely versatile since it is a general 
formulation of an inhomogeneous transmission line problem 
where metal conductors of any shape are placed upon one 


side of a dielectric substrate. 
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APPENDIX A 


Derivation of Spectral Domain Dyadic Green's 
Function Components (Open Boundary) 


Equations (2-3) and (2-5) are substituted into the 
continuity conditions (2-14) and (2-15) and the result is 


transformed into the Fourier domain which yields 


at y = 0 
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where иу =н, = H3 =u, and Е] = Е. = En so that 


The solutions to the two Helmholtz equations given by 


eguations (2-12) and (2--13) are now substituted; the 


h 


dependence of the coefficients As Еа Он under- 


stood and from here on is not shown explicitly. 


u A 2 е = 
c2 С“ = К со D (A-3) 
E - juuy B^ " -jayD" - jouy4D? (А-4) 
РИЕТИ он u 
Ko) C = co D (A-5) 
. h a h ; е u 
re + Jwe,y,B = -JjayD + Ju£g Y4D (А-6) 
р е æ 2_e_. 2-е = 
К со А“ = KB sinh y,D + К 26 cosh YD (A-7) 
-jayA” + шү АР = -jayB“ sinh Y4D = jayc“ cosh YoD 
и h ho 
-jouy4B cosh Y,D = Jwuy,C sinh YoD 
(A-8) 
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cO 2 
, е Е Du > 
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Dus DNE 2 Ame А = _ 
Kor K 458 Sinh Y,D КС cosh YD = J (a) (A-11) 


уд = jwe ny А“ + в sinh Y4D t jaye" cosh YoD 


DS cosh Y,D - c sinh YoD = J, (a) (A-12) 


The objective at this point is to replace the unknown 
coefficients AC through p? to obtain oa direct relation еве сеп 
the current density and electric field components. 


From equation (A-3) 


ре = (Es) * се 
cO 
and (A-5) 
k 
ph - (22) “ ch 
cO 


Substitute these in (A-4) and (A-6). Define the coefficients 





Such that 
ay >» 2 
F = went, — ЈЕ A-13a 
am UE ] ) 
Y k 
MEE =)“ (А-13Ь) 
12 сЗ 
then В? = FC? 4 са (А-14а) 
е = е һ 
апа Deere ЕС (A-14b) 
E, 2 Е2 1 
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Substitute (A-14) in (A-7) and define 


K Е 


с. = Е СИЕ sinh Y,D + cosh y,D) (A-15a) 
1 k "ND, 2 2 
cO 2 
Сы = = (62, 2 ы P D (A-15b) 
2 M mod Y2 
Ehen дз = сусе + G ch (A-16a) 


Likewise (A-8) can be rewritten if the following definition 
is used. 
Е 


IMMO. EE (D. \ = 
(С, E, Е sinh YD cosh Y4D) 


E e Fy cosh YoD (А= 156) 


= оа NS j 
G (6, + Би Fi sinh Y 4D) 


4 шиүі 
2 
= Y. (F5 cosh YD + sinh YD) (А-15а) 
then 
АП = бас“ + char (A-16b) 


Now solve for BS through ch in terms of A” and a 


Inversion of (A-16) yields 


С = ЈЈА“ + un (A-17a) 


С’ = JAS + ДА! (A-17b) 


181 





where 


> 4 _ 2 
а асс с 35 det С 
E det G 4 det G 
with det G = 616, - 6,65 


Substitute (A-17) in (A-14) and define 


0 u 
= — К Ј = ЕО 
1 Є, 2 1 E) 123 
€ 
0 u 
| и >к. 
2 E) 2,2 E, dl 
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e PJT Fag 


then 


DS KA" (A-19a) 


UJ 
| 


= К.А + KAP (A-19b) 


The next substitution is (А-17) апа (А-19) іп (А-11) and 


~ 


8-12). 
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Define 





Li = -(K, sınh YD + J, cosh YD) 
L, = 1 - E os inh D+ J h 
2 Koo 4 Sinh y, 4 COS YD) 
= оу 
L4 A (к. sinh YoD + J, cosh Y,D) 
a у 
HC. (к, cosh Y4D + Ji Sin Y4D) - l 
JE E 
E Su sinn y,Db 4 J, cosh YjD - 1 
4 €^ Y^ И о e d Тор ) 
- K, cosh YoD = Ј, sinh Y4D 
then 
2 e КЕ Be. _ 
K 201^ 4 К.012А = J (a) (A-20a) 
Hos. Y.L me ПЕ № АР = (10) (A-20b) 
013 DE DA 2 


Equations (A-20) may be solved for A 


h 4 e 
= and A 1n terms or 


the current density components. The result substituted in 


(A-9) and (A-10) will eliminate all unknown coefficients and 


express the electric field components in terms of the current 


density components with the following definition. 
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DM = (—— D EF 
E LE 2 E 
м. = (-29Y L о L y l 
1 n 2 ` 4 k 2. 3° DM 
co co 20 


u k 
C, ыл ы (_©2)^ L d. (ADAN) 


‚< 
| 





L 
_ 4 
M3 DM 
2 
4 QE AY» DM 
Then 
My M, J = 
м, M, J, E, 
If the M-matrix is inverted another set of coupled 
equations is obtained in the form 
Ny N, Е: 2. 
Е (A-23) 
N, N, E, J, 
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where with 


14 2 3 
М, “М, 
М. = МА = 
1 det M 2 det M 
(A-24) 
N = йш М = = 
3 det M 4 det M 


Note that the use of hyperbolic functions applies only 


if T О. вог Ze < 0 the trigonometric solutions for 


: ER Hnc 2 2792 
Ф, апа Y), apply since ү, = JY) Tor = К 2 о, 


155 





APPENDIX B 


Time Average Power Flow in the Spectral Domain 
(Open Boundary) 


The following expressions for the time average power flow 
in the three spatial regions are derived from equation (2-26). 
The solutions to Helmholtz's equations given in (2-12) and 
(2-13) are substituted in the transformed scalar potentials 
of (2-27). Although the detailed derivation would be quite 
lengthy and therefore prohibits a presentation here, the 
algebraic manipulation and the analytic integration with 
respect to y is routine. In the following derivation the 
dependence of the coefficients Af through JÉ on the variable 
a is understood. Note also that Үт “= Үз since regions 1 and 


3 are assumed air filled. 
In region 1 


Е а2+ү,2 
Р = =- — Ке Л BE u 


c 
un 


a 


* j2e(kg^A9*aAh - g/n?Ah*)] àa (B-1) 
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im region 2 for ү)? > 0, y, real 


1 E ug а+ү Р 
Р = - - Re = {sinh 2ү,р[ 


[2 
avg 4" 


ES Pe) 


2Y > 
ЕЕ 


2 
+ р(а?-ү,2) [е,(0[с®|2-|в® |2) + у([е? [2-8% 7 
N, е ех ех е 
* (cosh 2Y4D-1)———(€e,(B'C  -B C?) 
2 2Y, 2 
т рс Te у; 


+ 3 7 sinh o e 


E u Кы смс сд 


+ (cosh 2ү,р - 1) 182 (вёв\*+сёсћ*) 


-k* (88 *ph4c?*ch) 1) da 


(B-2) 
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in region 2 for 72 SUO oz 3Ү> ; Y> real 


р Da 
= = 1 WB : " i e DE e,2 
BE = – па ве a {== sin A ав ео) 


+ e 


+ D (a^*y27) te) (|B? | ^» |c? |) + не а ош ои 


козо h_hx* 


+ 2(1- сов 2730) go le, Im(B C^) * y Im(B'c^*)]) 
2 


hx_e hx 


san 2үзр{8° (В eos ce ee = 


) - k," (B cP-phc 


+ 
кој о 


+ 3 5 (1-cos 2үзр) {8°(в°в!*-с®еһ*®у „к? (в®*в!с®*сһВ, у) da 
(В-3) 
where Im( ) + "imaginary part of ( )" 
In region 3 
2 2 
со TY 
ЕТ ul 20 О 
Р ауд = gr Pe D [wB EE (ерір | ay ulD | ) 
* j2a(8^p?ph* - k,^p?*ph)] aa (B-4) 


The minus sign in front of all these integrals appears since 


a wave traveling into the -z-direction was assumed. 
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Bee slot transmission lines 


Observe that the coefficients BS Enrough p° are expressed 
in terms of A” and ge with equations (A-19), (A-17), (A-5) 
and (A-3). 

But A* and д? are expressed with (A-9) and (A-10) in 
terms of the transforms of the electric field components and 


are therefore known. 


For coplanar strip transmission line 


As for the slot line the coefficients B ehr och p? are 
expressed in terms of A* and rod (See above.) 

Since for this line the current densities are approximated, 
an E uicit dependence of these two ccefficients on the 
transform of the current density basis functions is required. 


This means equations (A-20) in Appendix A have to be 


inverted which yields 


M M 


e 3 ~ n 
А = => J, (a) + = J, fa) (B-5) 
k k 
El с 
Esel; k 3 L 
3 : 

Е о о в 

Kol EY DM WEY) Kol DM 


where all coefficients are as given in Appendix A. 
With this modification the expressions for the time 
average power (B-1) to (B-4) apply for the coplanar strip 


ime as well as the slot line. 


~ 
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APPENDIX C 


Perturbation Theory for Open Boundary Waveguides 


In this appendix a perturbational expression for the 
propagation constant of a wave traveling on an Open boundary 
waveguide structure is obtained. It is assumed the wave 
travels into the z direction and that the perturbation is a 
small change in the type of material within this structure. 

Let É and H represent the fields in the unperturbed 


Ү 2 


Structure where the z dependence as e is implied; then 


ҮхЕ- -jo[u] H (C-1) 
ў х Н = jwle] E (C-2) 
A 5 
and for the perturbed guide (z dependence e ја implied) 
V x É' = -jo[u*Au) É' (С--3) 
> а -> 
V x H' = ju[e+rde) E' (C-4) 


where "[ ]" denotes a tensor quantity. 
The last two equations are scalarly multiplied by H* and 
É* while the conjugates of the first two are scalarly 


multiplied by H' and É' respectively, which yields 
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DENEN o a? dt (C-5) 
El. Y x H* = -jufc]* Er - E' (C-6) 
H* . V x É' 2» -jo[y*Au] H' - H* (С-7) 
бк. ух Н! = jwletAc] É' - Еж (CSR) 


Equations (C-5) through (C-8) are now combined to obtain 


m. (8' x E*) 


Jwle+Ae] Е me ШӨЛ H* - H' (C-9) 


И 


ў. (Ёк х Ё!) 


-jw[e]* E* . Е' + ГЕЛИ] Н (C-10) 


where use of a vector identity is made: 


— 


namely ОКЕ хВ) 2-B- Vx A-A- Cx É 
Now add (C-9) and (C-10) and integrate over some volume V 
1 $7. Ш x È+) +7. (Их ВО] ау = 
ju n [[e*Ac] É' - É* - [cg] * E* . É'] àv 


-jw SJS Пи1% H* - H' - [ptapy) H' - É*] av (СОО 
V 


The restriction is now imposed that the material of the 


meaperturbed Structure iS isotropic so that 


ИЕСІ Е and similarly for [u] = [u]* 
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With this assumption some simplifications are possible. 
от ћег simplification of (C-11) is possible through an 


application of the divergence theorem. Then 


95 = 


Ji 


67 (8' х Ех + Не х Е!) 
S 


jw SSS [[A^c] * É' - E* + [Ay] Н' +» Н*] av (C-12) 
V | 
This equation is still exact and may be applied to any 


region V. Choose this volume such that a differential 


cylinder surrounds it, where the three sides are 


51 7 = а 
52 2 = Б = а+а? 
до 2 2 
53 wall of cyiinder at г = x Y 
Note that 
H' x Ех : n = -H' + n x É* 
and H* X Е! . n = -H* пх Е! 


For an open boundary structure and r very large the outward 


е > > 
unit normal vector n > r 


But lim |Н! х Вж + Н x Е! | ro 0 since the fields 


у -»ос 


finally vanish exponentially for large r. 
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This means that the surface integral over S, is equal to 


3 
zero. Denote the surface integration over Sy as F(a) and 


the surface integration over S, as (Б) then 
F(b) »-[F(a * $2 F(z)| 421. 
е а 


біпсе H' x EF* a«a НИ 5 


=> ! 
апа НЕ" се EE 


then 


Е. = a(y'ty*)] F(a) =; 


ОЕ 
< 2=а 


This derivation applies for any a along the z axis. 


Equation (C-12) can then be written as 


> 


у“) а> JJS (Н' х Ех + Н* x É') - n ds 
5 
= је /// [lade] E' Ё* + [Ay] H' . Не) ау (с-13) 
V 

where S is the cross-section of the structure and extends 
nouem finity. 

The right hand side of equation (C-13) remains essentially 
constant for a small enough variation in z so that the volume 
integral can be replaced by 


ju dz SS  [[Ae] É' * É* + [лди] H' - H*] as 
S 
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Poeing that [АЕ] = [et] [e] “ала [Ли] = [u'] - [Ij 


where [e] ES T КЕШЕ 


апа [u] = 99 (11) + (xD) 


the difference between the perturbed and unperturbed 


susceptibility tensors is found to be 


E po EE ee 


[^x] EOS» bu 


(C-13) can then be brought into the final form of 


-— 


ІШ / (EQ [Ax] Е ug [^x] . Н' . H*) ds 
Y аы = AS а ыс ш р, 
ff (B* x H’ 4+ Bt x H*) > ds 
S | Е 
The susceptibility tensors are internal tensors and 
therefore apply to the internal fields only. As it can be 
seen in Section 2.4.2 it is more convenient to use the 


Eterna! fields and external tensors if due to the physical shape 


of the ferrite demagnetization must be taken into account. 


=). __ © 
БЕДЕ ПЕЕ 
апа [x] H' = [у] В 


where the superscript "e" denotes external. 
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Equation (2-40) is now obtained. This detailed perturba- 
tional analysis of an open boundary waveguide could not be 
found in the literature and hence was rigorously carried 


out to prove its application for the slot line. 
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APPENDIX D 
Perturbational Expression for the Propagation Constant in the 


Spectral Domain 


Apply Parseval's relation to equation (2-42) which is 








then 
оно 2 aeos 27 
es = _—_—_————— ----- + (awe Ф 
A 915 9 Y, 
j MD eee 
- jawe, B (ð, $5 3m ЛИ 
2 оф 
T4 * их 
| + 2к-2Ху2 1528 3s Y3 + aue, 0,Y5] 
4 2 
VOLU A UA | ду da р) 


The transforms of the scalar potentials are known through 
equations (2-12) and (2-13) and the dependence of the 
coefficients BC ebreugh au on the assumed electric field 
distribution is derived in Appendix A. The integration with 
respect to y may be carried out analytically. Since this 
integration is performed in the spatial region 2, the two 
different solutions in this region for у" > 0 and w < 0 
have to be considered. The amount of algebraic manipulation 
is quite lengthy to obtain the final form of the integrand 
and will not be shown here, because it contains no new 


BEtormation. Note that in the lossless case the lefthand 


side of equation (D-1) is 
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ШЫ 559 bt = al + JAB 


A shorthand notation for the final form of equation (D-1) 


is then given by 


2 


+ eo р 2 


1 = " 2 
a = gn BD T | КА (B 


ауа -© 


Mr 2 


)? + (awe,)“ cc”) 


t2u€ oy, Bc? ] 


qe 2 h_h enh 
+ л C ł ошЕ С С] 


3 chy 24 


C23 о 


E Xz zh 1 


2 


+ 


Gc, EGv22 7 (CP) ^ €. (0862) ^ (89) ^. € 20e, ay, B*C^] 


2 Inm 4 má 2) 


eh " 
[*8y4B C + ашє,ВВ ] + у К (В) E 


4 


усу К 


DN 2 


и 2-h-h 2.0.8 
Ce С(ВҮу)) вс + (awe, ) BC 


+ 


+ we aby, (вев“ EEN 


ы 2 h, 2 lez eh hane 
O ESC О E owe (eB CC v Be )] 


4. 


+ k 4ghchj F | SIS 


X2z c2 


. (D-2) 


where the following notation is used: 


МЭЙ + 





The coefficients B through ch 


are either pure real or pure 
imaginary; if pure imaginary take the real part (see also 


Section 5) 


tor Y^ > 0 
5 = +1 
sinh 2 Y4D 
E, > + 
1 Ay, 2 
sinh 2 Y,D D 
E z 2 


cosh 2 Y4D - 1 


E == 
— —— 
3 Y) 


-- 


for Y^ О Үү = jys (replace Y> by Y> in (р-2)) 


4 = -l 

738 sin 2y5D , D 
1 45 2 

A EP sin 2Y5D 
Ba mn 


2 
1 - сос 2Y5D 
C3 — 7," — 
Y2 
furthermore, 


^ 


B. = Бе (у Ј 


198 








If the expression for AB” is desired, the following 


replacements have to be made: 


replace ү р ; 
р Хуу Y Хуу 


апа ya by X ' 


The superscript "+" on a' and A8 denotes that these 
parameters refer to a wave traveling in the positive z 
direction, so that the lower signs in the integrand have to 
be taken wherever a double sign appears; the superscript "-" 


refers similarly to a wave traveling in the negative 2 


direction where upper signs in the integrand have to be taken. 
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APPENDIX E 


Derivation of Spectral Domain Dyadic Green's Function 
Components (Closed Boundary) 


The continuity conditions given in equations (2-14) and 
(2-15) are transformed via (3-9) into the two dimensional 
Fourier domain. Transformed field quantities can be 
expressed through the solutions to the two Helmholtz's 
equations in equations (3-14)and (3-15) in forms of the 
coefficients Af through EM whose dependence on e and & 15 


understood from here on and is not shown explicitly. 


at y = 0 - 


| 


= Р h eV. : ME 
= 0, € - JwuY,B - 0, ED sinh YA) - Јоџмујр sınh У182 


(к=т) 


сё (2-52) - ре (k,^-£^) sinh Yıhz (E-2) 


c E h E _ D 
Jwe, Y,B OL, ES Jo£g Y D cosh YıR, а 5р cosh Y12) 


0564) 
ne) - D'(k,^-£7) cosh yh, (E-4 
at y =D 


e _. , Me: 
-a FA sınh У181 + jwuyjA sinh үл = 


CM e h 
E a, EB Sinh Үр - a,éC cosh YəP - Јоџув cosh YoD 


- jupy,C” sinh y4D (E-5) 
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е Da | е 26:2 e 2 2 
A (kg -E ) sinh Y4hi = ВО (К, -Е ) sinh y4D * C (Kk, iu ) cosh YoD 


(EEG) 
Е Е 
-a £A sinh YıRı + Jwyy,A sinh у181 + весе (E-7) 
e Ze 2 ' _ 
А (К) 56) sinh yh БИ Со (Е-8 


е ћ е 

~JWEQ YA cosh Y - Ay EA cosh YıRı - Jwe,y,B cosh YoD 
: біле. Ia h 

-2ше-ү-С sinh YD + a, EB sinh YD + a, EC cosh YoD 


= a) (E-9) 


А'(к„?-&°) cosh y¡h, - B"(k,2-£%) sinh y,D) 


1 
an 2152 
– С (К, -Е ) cosh үр = N) (E-10) 
Define 
e 
ee Sere о. ai 
(Ко -Е ) sinh у 185 
and 
ee 
2 = s ВЕБЕ o 
(Ko о соо У1 82 


then (E-2) and (E-4) can be rewritten as 


ре = ЕСД (E-12a) 


Jom h 
D = PoC 





The further objective is to express ве through ch in 


h which are related to the transforms of 


terms of Af and A 
the electric field and current density components at y - D. 
Substitute (E-12) into (E-1) and (E-3) and use the following 


definition: 


Сб, = a, & (Fy Sinh У1В2 - 1) 


Gy = Joy Y, F» sinh МЕР 


(E-13) 
G4 = JU£g YA F4 cosh МР 
Sas -ané (F3 cosh y,h, - I 2 
Then 
h Е h 

ЕЙ + СС” = jwuy,B 

| (E-14) 
ЕКЕ ссп = jwe y BÊ 

3 4 ve i 


(E-14) can be solved for a and ch in terms of БЕ апа gh 
to yield 
се = н ве + H.B” (E-16a) 
1 2 
в’ e h 
C = Н.В + H,B 
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with 


-Jwe,y,6, 
у 5 det G 

јоџу6, 
een 

Jo£,Y564 


Е 800 la) 


-jouY56G4 
ПИ астас 


with (E-16) now substituted into (E-5) and (E-6) and the 


@eeanition of 


E = -apé (sinh YD + H} cosh y,D) - jouy4H4 sinh ү,р 
J, = -jwuy, (cosh ФУ sinh YD) - a4 CH, cosh y,D 
(Е-1.7) 
J, = (К eo ey D + H, cosh D) 
3 2 Y2 il Y2 


28,2 
Ј, = Ну(к, -Е ) cosh y,D 


these are rewritten as 


e Dno e | ГТ 
71В + J4B = -А a8 sinh YıRı + JA WHY} Sinh Y1h4 


h 


e e 25.22 
J4B * J¿B A (ko -5') sinh У181 
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h 


where now ве and B can be solved for in terms of AS and 


AB. 
Define 
ПИ к у мул зз ои 
E n*4 Па 532 т 
1 det J 
M. Jouy, J; sinh МИ 
2 det J 
[a EJ, + (k 2 t?)g елп үп 
к = ne 0 lc т (Е-18) 
3 det J 
-jJupy 7J,sinh Ү1һі 
Ba = det J 
Then 
е _ e h 
ee A KA (E-19a) 
2 | 
в? = КЈА“ + ы : (E-19b) 
Now (E-16) and (E-19) are substituted into (E-9) and (E-10). 
Define 
L, 7 -jü€gY, cos yib > jwe,y,K] cosh y,D - joe, YA (Hi K, &H5K.) 


+ an K3 sinh YD + 0, 5 (H¿K, +H¿K3) cosh YD 
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ОЕ 22752 ови Ур WED Yo (HKD THK a) Sinks yop 


+ 0, EX y sinh үр + a, & (H3K,+H,K,) cosh YoD 


(E-20) 
ШЕ U-^-E4YK. sinh y.D - (k,^-E^)(H,K,4H,K,) cosh v4D 
3 zu З llo 2 ЗАВ Ү2 
Е 2 ae A DN | 
Ly = (ko -Е”) совҺ YP} (k, E IK, sinh y5D 
- {k CN КН КК.) сосһ y D 
2 зад 12 
Then 
о (Е-21а) 
1 2 2 
LAC+LAN= 7 (E-21b) 
3 4 > Ү? 
Equations (E-7) and (E-8) can be solved for АС апа Du 
yielding 
e Ez | 
um o0 ——— (E-22a) 
(Ko -£7) sinh Y1hi 
Е RE 
ЕШ - x n п 2 . (Е-22Ь) 


о. noon FM ОН T 
JL REM E ЕЯ juuy] (ko -E ) sinh Y3h4 


Substitute (E-22) in (E-21) and define 


Lio 
3 Зону sinn EN 
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L a. £L 
{ А n 2 


= 
|| 


4 2» | 2 
(ko -5") sinh Y1hi JWUY, (Ko -5") sinh hy 
Ly (E-22) 
M = 
аа хо 
1 JWHY, Sinn y,h, 
E E A+ Eg 
Е 
(Ko -E") sinh Yih4 juny, (ky -5") sinh Yihi 
Then 


ME + М.Е = J (E-23a) 


(E-23b) 


‚< 
Uy) 

m 
> 

+ 

= 
(24 

m 

|| 
Ly 
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ДРЕ ИРА F 


Time Average Power Flow in the Spectral Domain 
(Closed Boundary) 


A short hand notation for equation (3-34) is 


со 


Pavg = di БЕ Ке Pin E Pon Е P 3n! Ss 
where the first subscript denotes the spatiol region 1 to 3. 
These terms PE are then found by substitution of (3-17) 
and (3-18) into (3-34) where the variable £ is replaced by 
the constant 8. The integration with respect to y can be 
performed analytically since (3-17) and (3-18) show a simple 
dependence on this variable. Without showing the details 
of this rather lengthy but straightforward algebraic 
manipulation the single terms for equation (F-1) are given 
below with the following ее S Only a first order 
approximation will be used, hence Е, (а) = 0 and through 


(E-22a) A” (a) -0. 








Then 
2 2 
wunß сизү 
u 0 n 1 2__ 2 h 2 
a [ oe зу т (a Yi )h 1|A (о) | 


(852 
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2 2 








Pan = E jen zi sinh(2y4h4) – (а о -уј )ћ,) |ре (а ) |“ 
+ то" [ эу sinh (2y,h,) Е (a, ^-y, *)h5] | D^ (a, ) |^ 
+ а (kg^*8^) ре(а ) D'(o.) sinh(2y,h,) (F-3) 
For Pon it is not necessary to go through this cumbersome 


algebraic work if one observes that equations (3-17) and 
(3-18) are written according to the nomenclature used for 
the solutions for the open boundary Helmholtz's equations in 
the spectral domain given in equations (2-12) and (2-13). 
Therefore the derivations for region 2 in Appendix B hold if 
only the integrand is considered in equations (B-2) and (B-3) 
and a iS replaced by а. • 


eo for Y4^ О 


Р given by integrand of (B-2) dp 


2n 
where а, герігсев а 


ӘН Eor ү, < 0 


Ро, given by integrand of (B-3) (=) 
where un replaces a 


The coefficients An througl p^ for the above equations 
are related to the assumed transform of the electric field 


component Es as given in Appendix E with coefficients Еу to 
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K, and equation (E-22). All terms of (F-1) are then known 
and can be evaluated on the computer after a normalization 


technique has been applied which is described in Section 5. 
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APPENDIX G 


Normalized Equations for Open Boundary Lines 


In this appendix the equations of Appendix A leading to 
the dyadic Green's function components are rewritten such 
that quantities are scaled with respect to the thickness D 
of the dielectric substrate. The equations are first 
written out for the two cases where ү, > О апа Yo" < 0 
and are then combined in a form which lends itself towards 
a computer programming application. 

In the following notation an imaginary quantity u is 
written as ju" where u" is then real. The normalization 


— 


of the frequency with respect to D leads to 


D = l D 
ше во. 
€ 
Eun D 
“КОР ОХ 
«нор = 60 - (27) 2 


where now D/A is the normalized frequency and 


A 
XU 


то 
E 
|| 
№ 
24 
IU 
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From equations (2-1) and (2-2) then 


DS D 2 
(кау4б-- (2, 24 1094-11 
DE в. 2 E АХ ‚2 
(k 7D) = (27 7? les (т) | 
m 2 2 ; | А 
YP = [(aD)" - (К 10) ] which is a real quantity 


for ob Since l < ~ < EM 


pee case 1, ү) > 0 


According to equation (2-11) the above condition is 
D, 2 


о 
X (xr) ] so that 


equivalent to 00) ап [Е 


I 


za 


YD = lan) - ОО 15 гєа1. 


The coefficients Fy through Ny from Appendix A are then 











EN» 7) 
EE eb (0252 |, 
1 шыр у Kk ¿3D 
E n (Sc2?, 2 
2 YoD кур 
k_-D F 
Б. У: 2 ү 
E (x Б) = sinh y4D + cosh YD) 
el r 
E 9 


ZAG 





вр 5 2 
2201207) | 
BELL NE eee BY sinh ho 
2 к/р Бы d 2 
EU е, 
| F Ya D 
u" _ _aD BD г А Вр с и 
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M, = j M! 


3 3 
MY = LY/DM 
M = J M% 

2 
276 
det M = му му – му ми 
n; = Ue 
М = M¿/det M 
В м, 
N5 = -M5/det M 
Em м3 
М" = -M'/det M 
DES м" 
N^ E My /det M 
2 


b) Case 2, и о 


ЕО implies that the normalized transform 


variable (ар) < (21 2) 4 [Е == (5) 2] so that 


= 


У 


Y? = [ (k р) ^ - (aD)*]* 


214 





The coefficients Р. through N, are then 


K-D 
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DM = j DM" 
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k 1D 144 aD r 
М, - 7 Му 
М" ор BD al ne L ] zu 
l (кур) 4 E, YD Ka? 3° DM 
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ИО вр , (120m)? 17 (5625.2 ји] 1 
2 WED Y5D 2 ст Y5D k 1D l- DM 
0 | 
= L¿/DM" 
E На 
(К үр) 7 L 
M" = С 2 
4 o£ 4DY5D DM" m 


The equations for Ny through N, are identical with those 
listed above for a > 0 and need not be reiterated nere. 
For programming purposes it is now important to notice 
that indeed a complex arithmetic on the computer can be 
avoided since only a real arithmetic is necessary for the 
evaluation of the dyadic Green's function components M4 


through М4. А further simplification is even possible by 
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comparing the two sets of equations for the two cases 
I. > 0 апа е < 0 where it can be seen that some of 
tbe equations differ only by the hyperbolic or trigonometric 
sine and cosine function and some algebraic signs; the double 
prime superscript is immaterial for the further computation 
since it is only used to distinguish between real and 
imaginary quantities in the above derivation. It is then 
possible to reduce the two sets of equations above to one 
set where depending on whether s is greater or less than 
zero the proper signs and functions are applied. 

The procedure is self evident and needs no further 
development in this appendix. The result can easily be seen 
in the subroutine named "COMBO" of the attached computer 


program "COPLAN" where a consistent notation with the one 


used above is found. 
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АРРЕМОТХ Н 


Normalized Power Flow Equations 


The equations for the time average power flow of open 
boundary coplanar transmission lines are given in Appendix B. 
To accomplish a formulation which allows a numerical evalua- 
tion on the digital computer the following steps are to be 
taken. 

The coefficients AC through ph depend on the assumed 
field distribution at y = D as stated in Appendix А. These 
dependences are given for the two cases where 75 is greater 


or less than zero as follows in normalized form for slot lines: 


Yo“ > 0 

д^ В АП" 
p? 52 

и Е 
д^ P E, (0D) 
D ыр Үүр 
е п" 
sy = KR Ay 
D D 
ph _ pr" 
р2 р2 
ph" qn" 
= = К Zu 
D E D 
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9 Ee Yo. I = 
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xA oe 


| 

| 
q 
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; 





for Yo" < 0 


as above except 


р 
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With the assumption Е, (х) x 0 it follows that АЁ = 0, 


The normalized equations for the power flow then become 


Bor region 1 


2 2 
© (aD)~ + (y,D) h" 
Рр = - ат | gm жылт 297 n) GEM 


For region 2, Yo > 0 


E 1 TU 8D imi 2Y4D Ии °+(y,D)° 
p dm ED EINE ad ig 
EZ 
e e h" ch" 
+ (we pD( (5) * + e$ + OD . > E 
m | ЕСІ € p( (£5? = (2 E E m р( (2 52. Br 2 ) ] 
= “Ү2 ar Т H БЕСІ Е 
(ар) + (y,D) е се ве" с h" 
+ (cosh 2y,D - 1)——————— (w€£4D — + wuD —5- ) 
2 Y,D 2^ y 52 ЕС “т 
sinh 2y,D е Аһ" e ah" 
+ орг(вр) + (k9D) ^] [C— (+ Y) 
D D DW D 
cosh 2y-D-1 cout е АҺ" 
ПА | а (ар) (H-2a) 
D D D D 
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For region 2, йде < 0 





1 с2 8D sin 2y5D (ар) 2- (үлр) 
O T me" Dr D 
БЕ. D) 12 
el h h 
eS Er 22 Eu» 
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e gu e h" 
[ue D E)" + oo “| + өр BID) 2+ (кор) 7] 7 p | d (aD) 
(H-3) 





All coefficients in these integrands are real quantities 
now and can be easily programmed into a digital computer 
language. 

To accomplish an evaluation of the time average power 
flow for the coplanar strip line with real coefficients, 
equations (B-5) and (B-6) in Appendix B have to be written 
explicitly for the two cases Von > 0 and 1 < 0. 

Note that in the equations given here j, G0 20 is 
assumed. Furthermore, from equations (2-33) to (2-35) it 


is seen that the transform of the assumed current density 


J, (a) is imaginary so that 
J, (o) = j Jz (a) В 


Following the guidelines above and those in Appendix G, 


h 


the coefficients A* to р аге then derived in terms of the 


assumed current density. 


2 
For Yo > 0 
AS = M, T 
D^ (kp)? 2 
Сс 
an = 4 ADU 
p^ р2 
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e Mi 

= - J" 
p^ (kn)? Z 
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jh 7 Р А." 
p? p? 
E. қаш 
p? ü€4D Y5D K 7D DM" 
ве _ Р вё" 
р2 р2 
вё" _ T де А ^ ne 
p 1:2 2,2 
вр ae де m: дА" 
p^ 3 "i 4 p? 


and the rest of the coefficients are as given above for 


у" о. 


The normalized equations for the power flow are then, 


in region 1 


=> 2 
Em 1 BD (aD) + (y,D) e ne 
yn ae Le ST aco thy? + aww ey 
e h" 
+ aD [ (BD)? a (кор) 71 25 И | а (ар) H-4) 
D D 
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For regions 2 aná 3 equations (H-2) aná (H-3) apply 
without any further change so that here again a formulation 
with real quantities only Can be achieved. 

The further translation into a computer language is then 


routine. 
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APPENDIX J 


Normalized Equations for Closed Boundary Lines 


Sita lar, to the detailed manipulation of the various 
coefficients for the open boundary structures shown in 
Appendix G, a complex arithmetic for the closed boundary 
coplanar lines can be avoided. The equations leading to the 
transforms of the dyadic Green's function components are 
derived in Appendix E in a general form. It must be observed 
that the smallest value the discrete transform variable о 
can attain is zero. For the computation of the dispersion 
characteristic and subsequently for the line average power 
flow the transform variable t is replaced by the propagation 
constant B as stated in Section 3.2 where ki SS Ko Е 
so that for computations concerning these two line parameters 
Y, will always be positive but = may be positive or negative 
depending on the value of о, · 
In the computation for the resonant slot length, however, 
the transform variable E extends such that -œ < Е < о, 


thus in general three possible solutions for the transformed 


Helmholtz's equations (3-15), (3-16) must be considered, 


namely 
a) үу“ > 0 Р ee > 0 
b)- Y1^ O ; ла <0 
с) Y, < 0 , Y4* « 0 
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where a) and b) apply for the dispersion and characteristic 
impedance calculations as well. 

The coefficients in Appendix E for these three possible 
sets of solutions can then be derived in normalized form, 
where an imaginary quantity is written as u- ju". It is 
found that the transform of the dyadic Green's function 
component matrix [M] is imaginary and that the evaluation of 
each single component can be based upon an evaluation of real 
quantities only for the three cases mentioned above. 

Similar to the derivation indicated in Appendix G it is 
now possible to reduce these three sets of equations to one 
Ewhere the following notation is used for the three types 


of solutions: 


d 


Case a) upper signs apply 


Case b) middle signs apply, 'sinh YD" should be replaced by 
'sin ү5р' and 'cosh Y4D' by 'cos Y5D' ‚ furthermore 


І 
"ҮІ! replaces 'Y,D 


Case c) lower signs apply, replacements as indicated in 
b) above. Furthermore replace 'sinh y,h' by 'sin yyh' ; 


' cosh yP ' by 'cos yjh' and BD by 'yyD' 


The coefficients are then 


(кр)? ~ (р)? 
2_ 





+ 

mt 5 
<= = yy 1 1 
(Кр) (ED) ~] sinh МЕР 


228 





(k D)? ~ (8р)? 


F = I ERE 
2 Du 12 
((кор) -(£D2^] cosh ујћ, 
+ 
Су = aD ED e Pz Sinh МЕР - 1) 
әр 
с. = i WUD үүр P sinh Y185 
+ 
G, = М шє (р Уз О F4 cosh МЕР 
С, = - 4р 5р (Р. cosh Yıh, - 1) 
det G = G6, т G,6, 
H} = we,D y5D G,/det G 
+ 
+ 
Hz = E we,D Y,D G,/det G 
H, =  wuD Y,D G4/det G 
Jy - - aD ED (sinh үр + Н; cosh YD) + шур Y4D H4 sinh y,D 
= + 
J, = р шыр y,D (cosh YoD с Hy sinh YD) - o4D £D H, cosh y,D 
Ja = [(k D ND | (sinn Но в) 
3 72 2 1 2 
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2 2 
J4 = Н, [(k,D) = ер cosh YoD 
дек 7 - 2,7, - 2,2, 
EC В 2 2 P 
в = {+ a,D ED J} i [(XgD) ~ (5р) 13} Sinh уп) /аее J 
+ 


K, = у wuD үр Ja Sinh y¡h,/det J 
E 20 2 ! 
К. = H SADED ME [(k,D) (ED) 13, } Sinh yh, /det J 


Be m j WUD үүр 7; sinh y¡h,/det J 


Ly = к 6200 Үүр cosh МИ И ше-р YD Ky cosh Y,D 
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- шє,р YD (НК, - Н.К.) sınh YD E o, D ED K3 Sinh Y,D 


| + + 


+ a D ED (HK, HK.) cosh YD 


E Ер cosh У1 81 + шє,р ү,р К, cosh YD 


+ + 


Kk a, D ED (HK, - H,K,) cosh yoD 
L, = E [ (k De koo sinh yo D 
3 + 2 3 2 
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2 ДЕ ни | " 
s t H,K,) cosh yD 
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№ = [(к„р) 2- (ЕР) 2] cosh ујву + [(k4D)^- (ED) ^ JK, sinh yjD 


+ 


2 2 
[(k4D) - (£D) "] (H¿K, - H¿K¿) cosh y,D 


= 
|| 
t 


1 Е L,/wuD Үүр sinh y,h, 


+ 
— O nd 2 le 
I = | [L4 a D ED L,/wuD Y3D1/ [ (kgD) (ED) ^] sinh Y My 
M3 = Е L,/wyD Үүр sinh YM 
| 2 2 
M, = + [L4 - aD ED L,/wuD Y,D] ГА [ (kgD) -(ED) ] sinh МЕН 


While the above equations are applicable for the dispersion 
characteristic calculation of shielded slot lines and slot line 
resonator lengths, some further manipulations are necessary 
to provide means for a numerical evaluation of the tire | 
average power flow computation for slot lines. 


h 


The dependence of the coefficients nz through D оп the 


assumed electric field distribution is given in Appendix E. 
Rewriting these equations in the normalized form and combining 


the two sets of equations for the two possible solutions as 


2 2 


« 0 then leads to the following 


above where Y, > 0 ana Y» 


expressions. 


AT 
p* 


= 0 since in the first order approximation it is assumed 


ЕЕ 26 
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АВ E 

=. х 

й a 7 

D шур Үүр sinh үүр h7/D 
кА: 

p^? oUm p2 

me, ah 

р? - 4 р? 

+ ы Ве : ph 

р2 > 2 p? 

Е - в... > 

D^ 3 M 4 p? 

ре E се 
Ёз > 

mM ch 

c ES 


where the upper sign applies for Y4* > 0 and the lower sign 
For D e (Note that a quite similar notation was used 
for the expressions in Appendix H for the open boundary 
transmission lines, written SIS separately for и > 0 
апа < 0 while the above form shows the advantage of a 
more composite form.) 

The final step for a preparation for numerical analysis 


is then the normalization of the expressions for the discrete 


Fourier transformed power density terms given in Appendix F. 
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For region 1, 


h Sinh 2y,Dh,/D 
Е" БЮ А, 2 2 2 al 
Pin = wyuD КО E [Los D + (y,D) ] EET 
h 
T E LoD) ’- lD’) 


си тог region 3, 


L вр pe 2 | sinh 2y,Dh./D 
Pan - шедр > s | (ар + (yD) ] QR MM o 


22 72ү1Б 
Š T2 MapD) -vaD 23} 
+ оир Эр Ей С ^. (y p) 73 E 
kd Lo 27 (y, 021] 
Hs [ (BD) “+ (k gD) *] 2: >, sinh 2y,Dh,/D 


The two expressions for region 2 for oe > 0 and ү, < 0 
are readily available with the integrands of equations (H-2a) 
and (H-2b) in Appendix H due to the same nomenclature used 
for the solutions to the Helmholtz equations for open and 
closed boundary slot lines in equations (2-12, 2-13) and 
(3-17,3-18) where the continuous variable ap is replaced by 


the discrete variable aD. 
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The necessary translation into a digital computer language 
is then routine and can be easily seen in the subroutine "COMBC" 
of the attached computer program "RESO" which computes 
parameters for shielded slot lines and resonators. А 


consistent nomenclature with the one used above has been used. 
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APPENDIX K 


Circular Resonating Aperture on a Dielectric Substrate 


This appendix presents an attempt to analyze a circular 
slot on a dielectric substrate enclosed with conducting 
cylindrical walls of a circular waveguide closed at both 
ends by conductors across the cross-sections of this guide. 
A schematic of this configuration is shown in Figure K/l. 

This geometry is considered to be analyzed if depending 
on the assumed field mode of operation for a certain set of 
geometric parameters the resonant frequency can be found. 

The following derivation is not meant to be complete and 
rigorous but is merely intended to outline possible problems 
an investigator may encounter. For analytical purposes the 
conductors and the dielectric are assumed to be ideal and 
the dimensions of the surrounding walls are such that a 
regular circular cavity mode cannot exist. 


It is a well known fact that the scalar potential functions 


satisfy the Helmholtz equation which is in cylindrical 


coordinates 
19 (,39 9, 1 22% Г ao DENN (К-1) 
NES 0 225052 92 | 


If a general formulation of the resonating problem via a 
transformation is sought — following the guidelines in 


Section 2 and 3 — the partial differential equation (K-1) 
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Figure К/1 Circular Resonant Aperture on a Dielectric 
in a Closed Cylinder 
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is to be transformed into an ordinary differential equation. 
This suggests the use of a two dimensional transform which 
consists of a finite Fourier transform in the variable 0 and 
a finite Hankel transformation in the variable r. 


Thus define 


ШО 
no 
Н(п,5 ,2) = J eu re B (£r) h(0,r,z) dr de (K2) 
where В, is the Bessel function of the first kind of order n. 
It is easily shown that this transform applied to 

equation (K-1) will transform the differential operators of 
the variables 0 and r into algebraic operations so that an 
ordinary differential equation in the variable z remains. 

2 
Obviously the differential operator = is replaced by -n* 

од 

so that equation (K-1) via (K-2) becomes 
г EB. (E reise (r 22, LX чари на пи ан даг = 0 (K-3) 
0 по от 2 5а or 22 222 
where h is now a function of n, r and z. 


Integration by parts of the first two terms will then 


lead to 
a ов CTS 
oh me "m 
ЕЕ) х — gna en 
пег or О ox 0 
а oB_(E_r) a 
9 в м 2 1 ы 
Ша 55-5 / рза һаг 
a a 
(—- + kf) J r B (Ep) har = 0 (K-4) 
92. 0 
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Bessel's equation is known as 


У | Petia’ а 
5r t or =! 


2 
n 2 
YT Sm r) Ben) 


so that equation (K-4) can be written as 


„2 
СЕЛ Н(л,5,,2) = (Ё 


2 


2 
кон С) (К-5) 


which is the desired ordinary differential equation. 

However, this only holds if the boundary terms in (K-4) 
vanish. These terms must now be investigated in more detail. 
But first a sidestep is to be taken to formulate the 
behavior of the fielä components and their dependence on the 
scalar potential functions. The fields are contained within 

the circular waveguide cavity so that no energy is lost due 
to radiation. This fact makes the problem symmetricai in the 
9, r and z coordinates because standing waves do exist in 
these three directions. But then the fields may be expressed 
as TE or TM to any one of these coordinates where a conven- 
tional choice suggests the use of the z coordinate, similar 
eemevlindsical cavities. 

According to Collin [Ref. 14] one may assume a magnetic 


Hertzian potential UN = а, v(0,r,z) and an electric 


Hertzian potential Te a о, mete hull vor hybard 


2 


mode analysis. 
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The field quantities in the three cylindrical coordinates 


can then be expanded by Maxwell's equations to yield 


2 
E | 1 оф 
Br 7 rag ^ 200 г ур 
2 
IO . 9U 
Ба Ре та > ЈОН јр 
2 
Же 2 
Ea | E IO 
92 
(К-6) 
2 
E ж Ж 
Ho = groz ег ур 
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an 
|| 


К 5 55927 275 се 


2 
9 2 

С => 
02 


d 
И 


Bessel functions of the first kind are finite at the 
origin and have finite derivatives; similarly the scalar 
potential Bene for this eonfiguratıen ~ 1! be [Inite 
with finite derivatives so that the boundary terms in equation 
(K-4) vanish at the lower limit. At the upper limit r=a 
where Е, апа En must vanish (tangential electric field 
vanishes at the surface of a conductor), one of the boundary 
conditions for the scalar potential by equation (K-6) is 
then ф(а,0,2) = 0. This requires then that в (= а) = 0 
so that a discrete set of transform variables Em witch 
m= 1,2,3,... is obtained where a a is the m-th zero of 


the n-th order Bessel function, 
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while this derivation holds for the transformation of 
the potential function 6 it is not obvious that the same 
set of transform variables al apply for the transformation 
er the function y for which the boundary condition is ы = 0. 
A further complication in this analysis can be foreseen 
as follows. For the matching of the boundary conditions at 
the two interfaces at z = 0 and 2 = D the field components 
in the r and the 6 direction must be transformed where it is 
observed that the field components in the direction consists 
of partial derivatives of the scalar potential functions with 
respect to the variable r. The finite Hankel transform of 
a derivative of a function is given by Sneddon [Ref. 25]. 
For instance, if n ande E 


a 
Sp (E) = fx By (Ege) g(x) dr 


then 


|: г Ви (ЕЕ) ye аг 7 у=: SmOn-i 9g) 7 Uu бања (Em) 
It is seen that a derivative in the space domain is thus 
transformed into a algebraic operation in the transform 
domain but the order of the transforming Bessel function has 
to be changed which further complicates this analysis. 

As one may observe at this point, the relatively straight- 


forward transform method in a rectangular coordinate system 


is anything but simple in a problem of cylindrical geometry. 
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Finally, it should be mentioned De с Hertzian 
potentials are assumed to be in the other directions the set 
of equations (K-6) is changed in such a way that the boundary 
conditions are even more difficult to handle, 

This problem was abandoned here, since it was felt that 
if this method of analysis is applicable at all further study 


would more appropriately be the subject of new thesis work, 


241. 





APPENDIX L 


Design Curves 


This appendix includes some design curves for single 


slot and coplanar strip transmission lines. The appropriate 


parameters can be found in the figure captions. 
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